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ABSTRACT 

We present the first public data release (DRl) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. It consists of science- 
grade optical datacubes for the first 100 of eventually 600 nearby (0.005 < z < 0.03) galaxies, obtained with the integral-field 
spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto observatory. The galaxies in DRl already cover a wide 
range of properties in color-magnitude space, morphological type, stellar mass, and gas ionization conditions. This off'ers the potential 
to tackle a variety of open questions in galaxy evolution using spatially resolved spectroscopy. Two diff'erent spectral setups are avail- 
able for each galaxy, (i) a low-resolution V500 setup covering the nominal wavelength range 3745-7500A with a spectral resolution 
of 6. OA (FWHM), and (ii) a medium-resolution VI 200 setup covering the nominal wavelength range 3650-4840A with a spectral 
resolution of 2. 3 A (FWHM). We present the characteristics and data structure of the CALIFA datasets that should be taken into ac- 
count for scientific exploitation of the data, in particular the eff'ects of vignetting, bad pixels and spatially correlated noise. The data 
quality test for all 100 galaxies showed that we reach a median limiting continuum sensitivity of 1.0 x 10~^^ erg s~^ cm~^ A~^ arcsec"^ 
at 5635A and 2.2x 10~^^ erg s~^ cm~^ A~^ arcsec"^ at 4500A for the V500 and V1200 setup respectively, which corresponds to limiting 
r and g band surface brightnesses of 23.6 mag arcsec"^ and 23.4 mag arcsec"^, or an unresolved emission-line flux detection limit of 
roughly 1 x 10~^^ erg s~^ cm~^ arcsec"^ and 0.6 x 10~^^ erg s~^ cm~^ arcsec"^, respectively. The median spatial resolution is 3V7, and 
the absolute spectrophotometric calibration is better than 15% (Icr). We also d escribe the available interfaces and tools that allow 
easy access to this first pubhc CALIFA data at |http : //calif a . caha . es/DRl[ 
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1. Introduction 

^ The Calar Alto Legacy Integral Field Area (CALIFA) survey 
^ ( [Sanchez et al.|2012a| hereafter SI 2) is an ongoing large project 
of the Centro Astronomico Hispano-Aleman at the Calar Alto 
observatory to obtain spatially resolved spectra for 600 lo- 
cal (0.005 <z< 0.03) galaxies by means of integral field spec- 
troscopy (IFS). CALIFA observations started in Jun e 2010 with 
the P otsdam Multi Aperture Spectrograph (PMAS, |Roth et aL 



2005 ), mounted to the 3.5 m telescope, utilizing the large (74'' x 



64'') hexagonal field-of-view (FoV) oflTered by the PPak fiber 



* Based on observations collected at the Centro Astronomico His- 
pano Aleman (CAHA) at Calar Alto, operated jointly by the Max- 
Planck-Institut fiir Astronomic (MPIA) and the Instituto de Astrofisica 
de Andalucia (CSIC) 



bundle ( [Kelz et al.|[2006l [Verheijen eraLl[2004| ). A diameter- 
selected sample of 939 galaxies were drawn fr om the 7th data 
release of the Sloan Digital Sky Survey (SDSS, [Abazajian et al. 
20Q9| ) which will be described in Walcher et al. (in preparation), 



hereafter W12. From this mother sample the 600 target galaxies 
are randomly selected. 

Combining the techniques of imaging and spectroscopy 
through optical IFS provides a more comprehensive view of in- 
dividual galaxy properties than any traditional survey. CALIFA- 
like observations were collected during the feasibility studies 
( |Marmol-Queralt6 et alipOT l ; Viirone n"eFal.|2012|) and the PPak 
IFS Nearby Galaxy Survey (PINGS, Rosales'Ortega et al. |2010] ), 
a predecessor of this survey. First results based on those datasets 
already explored their information content (e.g 
[eraL]|2QT2l [Rosales-Ortega et al.|[20Tr 
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Fig. 1. Distribution of galaxies in the CALIFA mother sample on the 
sky (upper panel) and as a function of right ascension (lower panel). 
Galaxies in the CALIFA DRl sample are highlighted by filled symbols 
with red color. The number distribution in bins of 30° along the right 
ascension is shown in the lower panel for the mother sample (gray area) 
and the DRl sample (red shaded area). 



2011112012^1 ). CALIFA can therefore be expected to make a sub- 
stantial contribution to our understanding of galaxy evolution in 
various aspects including, (i) the relative importance and conse- 
quences of merging and secular processes, (ii) the evolution of 
galaxies across the color-magnitude diagram, (iii) the effects of 
the environment on galaxies, (iv) the AGN-host galaxy connec- 
tion, (v) the internal dynamical processes in galaxies, and (vi) 
the global and spatially resolved star formation history of vari- 
ous galaxy types. 

In this article, we introduce the first data release (DRl) of 
CALIFA which grants public access to high-quality data for a set 
of 100 galaxies. The properties of the galaxies in the DRl sam- 
ple are summarized in Sect.|2] We describe the data characteris- 
tics (Sect. [3]), data structure (Sect. [4]), and data quality (Sect. |5]) 
of the distributed CALIFA data as essential information for any 
scientific analysis. Several interfaces are available to access the 
CALIFA DRl data, which are introduced in Sect. [6] 

2. The CALIFA DR1 sample 

A sample of 939 potential CALIFA galaxies, also known as 
the "CALIFA mother sample", was drawn from the photomet- 
ric catalog SDSS DR7 as outlined in S12 and W12. The primary 
CALIFA selection criterion was the angular isophotal diameter 
(45'' < D25 < 80'') of the galaxies, which was complemented by 
limiting redshifts, 0.005 < z < 0.03. These redshift limits were 
imposed to ensure that all interesting spectral features can be ob- 
served with a fixed spectral setup of the instrument and that the 
sample would not be dominated by dwarf galaxies. Redshift in- 
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Fig. 2. Upper panel: Distribution of CALIFA galaxies in the w - z vs. 

color-magnitude diagram. Black dots indicate galaxies in the CAL- 
IFA mother sample (SI 2, W12) and colored symbols denote CALIFA 
DRl galaxies. Diff'erent colors represent the morphological classifica- 
tion of galaxies which range from ellipticals (E) to late-type spirals (Sd). 
Lower panel: The fraction of galaxies in the DRl sample with respect 
to the expected final CALIFA sample distribution in bins of 1 mag in 

and 0.75 mag mu - z. The total number of galaxies per bin in the 
DRl sample is written in each bin. The bins for which the number of 
galaxies in the mother sample is less than 5 are prone to low-number 
statistics and enclosed by a magenta square for clarity. 



formation was taken from SIMBAD for all galaxies where SDSS 
DR7 spectra were unavailable. The reader is referred to W12 for 
a detailed characterization of the CALIFA mother sample and a 
thorough evaluation of the selection efifects implied by the cho- 
sen selection criteria. From the CALIFA mother sample, 600 
galaxies are randomly selected for observation purely based on 
visibility, and we refer to these galaxies as the virtual final CAL- 
IFA sample hereafter. 

The first 100 public DRl galaxies were observed in both 
spectral setups from the start of observations in June 2010 un- 
til June 2012. We list these galaxies in Table [T] together with 
their primary characteristics. The distribution of galaxies in the 
sky follows the underlying SDSS footprint (Fig. |T]). An excep- 
tion is the additional cut in ^ > 7° that was applied to the region 
of the Northern Galactic Cap, given the sufl&ciently large number 
of objects with better visibility from Calar Alto. The number of 
galaxies in DRl is not homogeneous as a function of right ascen- 
sion, Qf(J2000), and has a clear peak around a ~ 255°. This can 
be explained by the unexpected downtime of the 3.5 m telescope 
from August 2010 until April 201 1 due to operational reasons at 
the observatory, which delayed the survey roughly by 8 months. 
CALIFA observations, therefore, so far span three summer sea- 
sons at Calar Alto, but only a single winter season. Neverthe- 
less, the distribution of physical properties is nearly random, as 
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Fig. 3. Luminosity functions in the r band of the DRl sample (red 
points), the CALIFA mother sample (blue points), and the low-redshif t 
sample of the NYU value added catalog (black line,| Blanton et al.|2005| ). 
Despite the small size of the DRl sample the luminosity function is 
reproduced well. 



expected, and covers galaxies with a wide range of properties as 
discussed below. 

The distribution of galaxies in the color-magnitude diagram 
(Fig. [2]) shows that the DRl sample covers almost homoge- 
neously the full range of the CALIFA mother sample. On aver- 
age, the DRl targets comprises -20% per color-magnitude bin 
of the total expected number when CALIFA is completed. How- 
ever, there is currently still a deficit of low luminosity galaxies 
with intermediate colors. In other color-magnitude bins, espe- 
cially within those where the mother sample contains few galax- 
ies, fluctuations can be explained by the efl'ect of low number 
statistics. Figure |2] highlights the need to further increase the 
numbers to the full CALIFA sample to obtain enough galaxies 
in each bin for a meaningful multi-dimensional statistical anal- 
ysis. In Fig. [3j we show the r band luminosity function (LF) 
of the DRl sample as compared to the mother sample and the 
reference SDSS sample of Bla nton et al.| ( |20Q5] ). All technical 
details on how we obtained the LFs are described in W12. We 
simply note here that despite the small sample size of DRl, we 
already reproduce the LF reasonably well. The turnover of the 
LF at > -19.5 is entirely expected and understood. 

Galaxy morphologies were inferred by combining the inde- 
pendent visual classifications of several collaboration members 
as described in W12. In Fig. |4] we show the fraction of DRl 
galaxies with respect to the expected final sample distribution for 
diff'erent morphological types grouped into elliptical, lenticular, 
spiral galaxies (separated by difl'erent bar strength) as well as on- 
going galaxy mergers. A more detailed classification of spirals 
into early- and late-types is available, but we do not distinguish 
between them here because of the still modest number of galax- 
ies within DRl. The DRl coverage seems to be consistent with 
a random selection because the fraction of DRl galaxies with 
respect to the expected final sample is almost constant for all 
types. Axis ratios (hid) were measured from the SDSS r band 
image by calculating light moments after proper sky subtraction 
and masking of foreground stars (see W12 for details). The frac- 
tions of axis ratios, which can be used as an indication of the 
inclination of spiral galaxies, covered by the DRl sample is ho- 
mogeneous with respect to the final sample (Fig. |5]). Based on 
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Fig. 4. The fraction of galaxies in the DRl sample with respect to the 
expected final CALIFA sample distribution, split by visually classified 
morphology. We divide the galaxies into ellipticals (E), lenticulars (SO), 
non-barred spirals (SA), weakly barred spirals (SAB), strongly barred 
spirals (SB), and ongoing mergers (M) of any type. The morphological 
distribution of the DRl sample lies close to that of the mother sample. 
The total number of galaxies in the DRl for each morphology type is 
written on the bar. Error bars are computed from the Poisson errors of 
the associated DRl number counts. 
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Fig. 5. The fraction of galaxies in the DRl sample with respect to 
the expected final CALIFA sample distribution, as a function of axis 
ratio {bla). Galaxies were separated into early-type galaxies (E+SO) 
and late-type galaxies (Sa and later). The CALIFA mother sample does 
not include any elliptical galaxies with b/a < 0.3 or any spiral galaxies 
with bla> 0.9. Error bars are computed from the Poisson errors of the 
associated DRl number counts. 



a Kolmogorov-Smirnov test, we further quantified that the mor- 
phology and the axis-ratio distribution of DRl is consistent with 
being randomly drawn from the CALIFA mother sample with 
>95% confidence. 

In Fig. [6j we present the distribution of stellar masses for 
the DRl galaxies. They were directly inferred from the CAL- 
IFA data by applying spatially resolved spectral synthesis us- 
ing the Starlight code (Cid Femandes et al.|2QQ5| ) with a com- 
bination of single stellar populations (SSPs) from the libraries 
of [Gonzalez Delgado et aL] ( [2005] ) and jVazdekis et al.[ ( [2QT0| ), 
which both adopt a Salpeter initial mass function. Details of 
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Fig. 6. Distribution of stellar masses for galaxies in the DRl sample. 
The stellar masses have been determined from the CALIFA data (see 
text for details). 
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Fig. 7. Emission-line diagnostic diagram for the most central spaxel of 
each DRl datacube. Only objects for which all re quired emission line s 
have S/N > 3 are show n. The demarcatio n lines of|Kewley et al.lpOOT] ), 
IKauff mann et al.| ( [200 3 ), and C id Fernandes et al.| |2010') are used to 
classify the galaxies into star forming (SF), intermediate, Seyfert, and 
LINER-type galaxies, which are denoted with black, magenta, red, and 
green symbols, respectively. 



these measurements will be presented in Gonzalez Delgado et al. 
(in preparation). The DRl galaxies cover intermediate to high- 
mass galaxies, including more than 5 galaxies per 0.25 dex bin 
between 10^^ and lO^^M© and a median value close to 10^ ^M©. 
The asymmetric distribution is expected from the distribution in 
absolute magnitudes (see Fig.|2]) and is inherited from the CAL- 
IFA mother sample due to its selection criteria (see W12 for de- 
tails). 

Many different kinds of galaxies are already covered in the 
DRl sample in sufficient numbers to perform spatially resolved 
comparison studies. In addition to global galaxy parameters 
presented previously, we show the [Oni]/I^ vs. [NnJ/Hof 
emission-line diag nostic diagram ( [Baldwin et al.|198l] |Veilleux 



|& Qsterbrock||198 7 ) for the nucleus in Fig. [T] It is constructed 
from the most central spaxel at the coordinates of each CALIFA 
galaxy center (Table [T] excluding objects with a wrong astrom- 
etry in the data as discussed later in Sect. |5.2| ). After removing 
the stellar continuum using a spectral synthesis approach, we 
measured the emission line fluxes by fitting Gaussians to their 
profiles in the residual spectra. A variety of diff'erent ioniza- 
tion mechanisms can be found in the DRl galaxy sample, rang- 
ing from pure star formation to Seyfert-type active galactic nu- 
clei (AGN), with a significant number of galaxies in between. 
Also, Low-Ionization Nuclear Emission-line Regions (LINERs) 
([Heckman|| 1980| ), which are predominantly hosted by bulge- 
dominated or elliptical galaxies, appear to be frequent. Robust 
classification of individual galaxies according to such a scheme 
can be complicated, i.e., in border-line cases or due to system- 
atic uncertainties in the emission-line measurements, as well as 
adopting other types of diagnostic diagrams (e.g. Cid Fe rnandes | 
|et al.|[2Q10| ), or a diff'erent classification scheme. This has been 
considered in the [Kehrig et al.| pQ12| ) study of the ionization 
source of the interstellar medium in two LINER-like early-type 
CALIFA galaxies. Rather than to suggest or impose a certain 
classification of individual galaxy nuclei for the DRl sample, we 
present one particular representation here only to demonstrate 
the diversity of the DRl galaxy sample. 



3. Data processing and error propagation 

The instrument characteristics and observing strategy of the 
CALIFA survey define the requirements for the data reduction 
scheme. These requirements are thoroughly described in S12 
and only briefly summarized here for completeness. The PPak 
fiber bundle of the PMAS instrument comprises 382 fibers, each 
with a diameter of 2.7'' projected on the sky. The primary FoV 
has a hexagonal shape of 74''x62'' in size, sampled by 331 fibers 
with a filling factor of -60%. A set of 36 fibers are dedicated to 
sample the sky background and are distributed in 6 smaller bun- 
dles along a circle of ~ 72'' radius around the FoV center. The 
remaining 15 fibers are coupled to the calibration unit of the in- 
strument. 

Each CALIFA galaxy is targeted twice with diff'erent spec- 
tral settings, a low-resolution (V500) setup covering the nominal 
wavelength range 3745-7500A at spectral resolution of 7? ~ 850, 
and a mid-resolution (VI 200) grism covering the nominal wave- 
length range 3400-4840A at a spectral resolution of 7? ~ 1650. 
The useful wavelength range is, however, reduced by internal 
vignetting within the spectrograph for several fibers to 4240- 
7140A and 3650-4620A in the worst case for the V500 and 
VI 200, respectively. Three dithered pointings are taken for each 
object to reach a filling factor of 100% across the entire FoV. 
The exposure time per pointing is fixed to 900 s for the V500 
and 1800 s for the VI 200. The latter is further split into 2 or 3 
individual exposures. 

We are continuously upgrading the CALIFA pipeline (see 
S 12 for details). The main improvements to the data reduction 
pipeline used to produce the DRl data are briefly mentioned in 
the next section followed by a detailed characterization of the 
propagated noise. 

3.1. Improvements on the CALIFA data reduction scheme 

At the first stages of the project, we employed the R3D data re- 
duction package ( [Sanchez [20061 as the basis to develop a fully 
automatic IFS data reduction pipeline dedicated to the CALIFA 
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Fig. 8. Example of a stellar population fit with starlight to the central 
V500 spectrum of NGC 5966, restricted to the unvignetted part. The 
observed, modeled and residual spectrum are shown together with the 
pipeline produced ±\cr error. Sky and galaxy emission line regions 
masked during the fitting are indicated by the gray shaded areas. 



data. A description of that pipeline and the individual reduc- 
tion steps was presented by S12. With the aim of improving 
the flexibility, portability, maintenance and capabilities of the 
pipeline, most of the data reduction tasks have been transformed 
to a Python-based architecture. 

The main driver for the pipeline improvements was to imple- 
ment the propagation of the Poisson plus read-out noise as well 
as bad pixels caused by cosmic ray hits, bad CCD columns, or 
the efl'ect of vignetting, from the raw data to the final CALIFA 
data product in a consistent way. This was not possible with 
the reduction scheme adopted by the previous pipeline version 
(VI. 2), because the numerous interpolation and resampling steps 
did not allow a reliable error propagation. Details of the most 
important improvements for the current pipeline version (VI. 3c) 
are provided in Appendix [A| 

The final output of the VI. 3c pipeline are cubes of inten- 
sity, errors and bad pixels. Details on the data format are de- 
scribed in Sect. |4] Below we evaluate specifically the accuracy 
of the pipeline-produced noise estimates and characterize eflfects 
of correlated noise that should be properly taken into account in 
any scientific analysis. 



3.2. Accuracy of the propagated noise 

The derivation of error spectra in IFS data is a complex task, for 
which there are no well established recipes. It is therefore rel- 
evant to verify the reliability of the pipeline errors (e^). An ap- 
proximate assessment of the error spectra was carried out with 
the aid of full spectral continuum fitting described by Cid Fer- 
nandes et al. (in preparation), and Gonzalez Delgado et al. (in 
preparation). The CALIFA spectra are decomposed in terms 
of combinations of SSP spectra fr om [Gonzalez Delgado et al.| 
([2005) and Vazdekis et al. ( 2010|) inc luding a dust extinction 
term following the Cardelli et al. ( |1989| ) law. An example of such 
a fit is shown for a central V500 spectrum in Fig.[8]for illustrative 
purposes. For a certain spectrum k it presents the observation 
(0;ik), the best-fit model (M^/.), the residuals (R;ik = O^ik - ^Ak), 
and the provided pipeline error spectrum (6^^^). 
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Fig. 9. Histogram of the observed (O^^) minus model (M^^) residu- 
als at each wavelength normalized by the formal flux error (e^^t)- The 
model spectra are obtained from starlight fits for the stellar continuum, 
excluding the main emission lines, as well as bad pixels. The histogram 
describes almost a Gaussian centered at with a dispersion of ~ 0.8. 
The blue line indicates a Gaussian distribution centered at with a dis- 
persion of 1, normalized to the total number of pixels for comparison. 
On the assumption that the spectral residuals are due to noise (i.e., ne- 
glecting model imperfections), this figure demonstrates that the pipeline 
error estimates are reliable. 



Neglecting systematic deviations of the model spectra with 
respect to the real data in some wavelength regions, one would 
expect the residuals R;ik to be typically of the order of 6^^^. In- 
deed, we find that this is the case. The histogram of U;ik = 
Rxkl^xk (Fig -[21 was obtained from over 10^ data points of nearly 
10^ spectra for the V500 spectra distributed within this DRl. 
Emission lines and faulty pixels were excluded from the statis- 
tics (shaded areas in Fig. [8]), as the spectral fits are only meant to 
reproduce the stellar component. The histogram follows nearly 
a Gaussian shape, centered at -0.04 with a standard deviation 
of 0.8. The expected standard deviation would be LO if the er- 
ror spectra would statistically agree with the noise in the resid- 
ual spectra and assuming that the stellar population model is a 
perfect fit. Selecting diff'erent wavelength ranges for this com- 
parison to avoid regions with systematic template mismatches 
mainly changes the centroid of the distribution whereas the stan- 
dard deviation is nearly constant. Thus, the propagated errors are 
overestimated by -20% in the mean. This was also verified for 
the VI 200 setup using the same method. The overall conclusion 
of this test is that the errors provided in DRl are robust, with a 
slight systematic overestimation of -20% that can be taken into 
account during any analysis if needed. 

3.3. Characterization of spatially correlated noise 

It is often necessary to spatially co-add spaxels in the datacubes 
to achieve a minimum S/N in the spectra for a specific applica- 
tion. For CALIFA we adopt an inverse-distance weighted image 
reconstruction scheme so far, which averages the flux among all 
fibers within 5'' for a given spaxel in the final datacube by as- 
suming a 2D Gaussian profile with a dispersion of V for the 
individual weighting factors (see S12 for details). Like many 
other image resampling schemes it introduces significant corre- 
lation between the spaxels in the final datacubes. This can be un- 
derstood from simple arguments. Each CALIFA dataset contains 
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Here, we characterize the effect of the correlated noise by 
determining the ratio of the "real" error (6reai), directly estimated 
from the residuals to the analytically propagated error (6bin) 
of binned spectra as a function of bin for a certain target S/N. 



100 



Fig. 10. Ratio of the real error (6reai) to the analytically propagated 
error (6bin) as a function of number of spaxels per bin for all the V500 
(upper panel) and VI 200 (lower panel) data of DRl at a target S/N 
of 20. Grey shades mark the Icr, 2<t and 3cr levels. The blue lines 
represent the best fit logarithmic function with or = 1.38 and a = 1.35, 
respectively. 



993 physically independent spectra from the fibers of the three 
dithered pointings, but the final datacube consists of more than 
4000 spectra at a sampling of x per spaxel. This automat- 
ically implies that the final spaxels cannot be completely inde- 
pendent from each other within the CALIFA datacubes through 
the complex correlation of signal and noise between neighboring 
spaxels. 

In the limit case of co-adding the spectra of the entire cube, 
the pipeline analytically calculates an error weighting factor for 
each pixel such that the formal error of the co-added spectrum is 
identical to the one obtained by co-adding the individual fiber 
spectra. Of course, this is an unrealistic case, because only 
spectra within a small zone are typically co-added to preserve 
some spatial resolution. A popular method for adaptive bin- 
ning is the Voronoi-binning scheme implemented for optical IFS 
data by Cappella ri & Copin] ( |2QQ3| ). It assumes that the spec- 
tra are completely independent of each other to compute the re- 
quired bin size for a given target S/N. Blindly adopting such a 
binning scheme to CALIFA datacubes will lead to incorrect re- 
sults, because the assumption that the spectra are independent 
is not valid. Either the bin sizes will be smaller than required 
to achieve the target S/N, or alternatively expressed, the error of 
the co-added spectra will be higher than formally expected given 
the error of individual spaxels. 



The results obtained for all DRl datasets are shown in Fig. 10 
for the two instrumental setups with a target S/N level of 20. 
The observed trends can be sufl&ciently described by a simple 
logarithmic functior[^ 

^real = ^bin[l +a^l0g^] , (1) 

with n the number of spaxels per bin. 

The values for the slope a range from 1.35 to 1.45 with a 
mean of a = 1.4, for target S/N values between 10 and 60. The 
fits to the data result in less accurate exponents for target S/N 
ratios of 10 and 60 given the poor sampling towards large- and 
small- size bins, but the overall shape is preserved. This means 
that a mainly depends on the number of spaxels per bin and not 
on the target S/N. 

4. CALIFA data format and characteristics 



The CALIFA data are stored and delivered as datacubes (three- 
dimensional data) in the standard binary FITS format and consist 
of four FITS Header/Data units (HDU). These datacubes repre- 
sent (1) the measured flux densities, corrected for Galactic ex- 
tinction as described in S12, in units of 10"^^ erg s"^ cm"^ A"^ 
(primary datacube), (2) associated errors, (3) error weighting 
factors, and (4) bad pixels flags (Table [2]). They allow users to 
properly take into account the characteristics of each dataset con- 
cerning bad pixels and noise for their specific analysis. The first 
two axes of the cubes correspond to the spatial dimension along 
right ascension and declination with a V x sampling. The 
third dimension represents the wavelength along a linear grid. 
The dimensions of each datacube {Na, Ns, and A^^), as well as the 
spectral sampling (J^) and constant resolution (^^) along the en- 
tire wavelength range, are summarized in Table [3] Note that the 
spatial resolution is worse than the actual seeing during the ob- 
servation because the large aperture of each fiber (2V7) strongly 
undersamples the point spread function. The final CALIFA spa- 
tial resolution is mainly set by the dither pattern and the image 
reconstruction scheme rather than the seeing. We evaluat e thi s 
as part of our quality control tests discussed below in Sect. |5^ 



4.1. Error datacubes 

The error datacubes in the first FITS extension correspond to the 
associated Icr noise level of each pixel as formally propagated by 
the pipeline. As validation of our eflforts to optimize the pipeline, 
we verified that the measured noise is systematically lower by 
-20% than the formal error vector for individual spaxels (see 
Sect. |3.2| ). In the case of bad pixels, we assigned an error value 
that is roughly ten orders of magnitude higher than the typical 
value for the considered dataset. Any analysis based on the x'^ 
statistics will implicitly take into account bad pixels if the error 
vector is considered. The correlation of noise becomes important 
only in cases where the CALIFA data need to be spatially binned 
as discussed in Sect. |3.3[ The second FITS extension reports the 
error scaling factor for each pixel in the limiting case that all 
valid spaxels of the cube would be co-added. Suggestions about 
how to deal with the correlated noise and proper usage of the 
error scaling factors were described above in detail. 



^ An independent estimate of this empirical relation, proposed by Cid 
Fernandes et al. (in preparation), is completely consistent within the 
estimated uncertainties. 
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Fig. 11. Fraction of valid pixels in each spectrum across the CALIFA FoV (left panels) and as a function of wavelength (right panel) for the 
V500 and VI 200 setup. The spaxels most severely affected by the vignetting in the blue and red part of the spectra lead to the ring like structure 
around the FoV center. The blue part of the VI 200 data does not show a vignetting effect here because we cut off the wavelength range in the final 
data already at 3 65 OA due to the very low sensitivity in the blue. Four bad CCD columns are visible that lead to a significantly reduced fraction of 
valid pixels at narrow wavelength regions. 



Table 2. CALIFA FITS file structure 



HDU 


Extension name 


Format 


Content 





Primary 


32-bit float 


flux density in units of 10"^^ erg s"^ cm"^ A"^ 


1 


ERROR 


32-bit float 


Icr error on the flux density 


2 


ERRWEIGHT 


32-bit float 


error weighting factor 


3 


BADPIX 


8 -bit integer 


bad pixel flags (l=bad, O=good) 



Table 3. Dimension and sampling of CALIFA datacubes 



Setup A^^^ Ns^ 



V500 
VI 200 



78 
78 



73 
73 



1877 
1701 



3749A 
3650A 



7501A 
4840A 



2.0A 
0.7A 



6.0A 
2.3A 



Notes. Number of pixels in each dimension. Wavelength of the first 
pixel on the wavelength direction. Wavelength of the last pixel on the 
wavelength direction. Wavelength sampling per pixel. Homoge- 
nized spectral resolution (FWHM) over the entire wavelength range. 



4.2. Bad pixel datacubes 

Bad pixel datacubes are produced by the pipeline and stored in 
the third FITS extension. They report pixels for which no suffi- 
cient information is available in the raw data because of cosmic 
rays, bad CCD columns, or the effect of vignetting. The data at 
bad pixels have been interpolated, and it is considered not usable 
even if the spectrum looks "good" at those locations. The VI 200 
data is less affected by cosmic rays compared to the V500 data, 
because several frames are observed per pointing. The data at a 
bad pixel are therefore restored from the other available frame(s) 
although this results in a lower S/N. 

The distribution of bad pixels is not homogeneous within 
the datacube because of the vignetting effect as noted in S12. 
In Fig. [TT] we present the typical fraction of valid pixels along 
the spatial and spectral axis for the V500 and VI 200 setup, re- 
spectively. Note that four bad CCD columns have been identi- 
fied, which lead to four wavelength regions where the fraction 
of valid pixels is significantly reduced. The vignetting effect of 
the instrument results in a wavelength coverage reduced by up to 
-25% for some fibers on the blue and/or red side of the spectral 
range. 



The hexagonal PPak FoV is resampled to a rectangular grid, 
-SO that the uncovered corners are filled with zeros. These are also 
flagged as bad pixels for consistency, whereas the residuals of 
bright night-sky emission lines are not flagged as bad pixels. The 
strength of their residuals is different for each dataset and might 
"t)e handled differently depending on the specific data analysis. 

4.3. FITS header information 

The FITS header contains the standard keywords that register 
the spatial axes to the standa rd World Coordinate System (WCS, 
Greisen & Calabretta||2002| ) and the wavelength to the spectral 
axis in a linear grid. Each CALIFA datacube contains the full 
FITS header information of all raw frames from which it was 
created, where each header entry is expanded with a unique pre- 
fix for a given pointing/frame. The prefix consists of the string 
"PPAK" followed by the designation of the pointing PREF, 
which can be "PI", "P2", or "P3" for the three pointings of the 
V500 setup, or "PlFl", "P1F2", "P2F1", etc., for the individual 
frames taken for each pointing of the VI 200 setup. 

The reduction pipeline also collects information regarding, 
e.g., sky brightness, flexure offsets. Galactic extinction, approx- 
imate limiting magnitude, etc., and adds it to the FITS header. 
Header keywords that may be of general interest for the anal- 
ysis and/or evaluation of the data are summarized in Table |4] 
for convenience. Note that the systemic velocity of the galax- 
ies (MED_VEL) automatically estimated by the pipeline is not 
robustly measured from the VI 200 data, because of the smaller 
wavelength coverage and lower S/N compared to the V500 data. 



5. Data Quality 

This first CALIFA data release provides science grade data of 
100 galaxies to the community. Hence, in order to define the 
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Table 4. Main FITS header keywords and their meaning 



Keyword 


data type 


Meaning 


OBJECT 


string 


Name of the target galaxy 


CALIFAID 


integer 


ID of the galaxy 


NAXISl 


integer 


Number of pixels along right ascension axis (A^^,) 


NAXIS2 


integer 


Number of pixels along declination axis (A^^) 


NAXIS3 


integer 


Number of pixels along wavelength axis (A^^) 


CRPIXl 


float 


Reference pixel of the galaxy center in right ascension 


CRVALl 


float 


Right ascension a (J2000) of the galaxy center m degrees 


CDELTl 


float 


Sampling along right ascension axis in arcsec 


CRPIX2 


float 


Reference pixel of the galaxy center in declination 


y^T^'¥ TAX 

CRVAL2 


float 


1* A ' c / ^ ^ r\r\r\\ r*ii i a • i 

Decimation o (J 2000) of the galaxy center in degrees 


CDELT2 


float 


Sampling along declination axis in arcsec 


CRVAL3 


float 


Wavelength of the first pixel along the wavelength axis in A 


CDELT3 


float 


Sampling of the wavelength axis in A 


hierarch PIPE VERS 


string 


Version of the reduction pipeline 


hierarch PIPE UNITS 


string 


Units of the flux density 


hierarch PIPE GALEXT AV 


float 


Galactic V band extinction along line-of-sight 


hierarch PIPE REDUDATE 


string 


Date/time the data were reduced 


hierarch PIPE VMAG 3SIG 


float 


Estimated 3cr surface brightness depth in mag arcsec"^ 


MED_VEL 


float 


Estimated systemic velocity in km/s 


hierarch PPAK PREF DATE-OB S 


string 


Date/time of the observation 


hierarch PPAK PREF UT_START 


float 


Start of the observation in seconds UT time 


hierarch PPAK PREF UT_END 


float 


End of the observation in seconds UT time 


hierarch PPAK PREF MJD-OBS 


float 


Modified Julian date of the observation 


hierarch PPAK PREF AIRMASS 


float 


Airmass during the observation 


hierarch PPAK PREF EXPTIM 


float 


Exposure time of the observation 


hierarch PPAK PREF HVCOR 


float 


Line-of-sight diff'erence to helocentric velocity 


hierarch PPAK PREF EXT_V 


float 


V band atmospheric extinction in magnitudes 


hierarch PPAK PREF PIPE FLEX XOFF 


float 


Measured flexure ofl'set in x-direction in CCD pixels 


hierarch PPAK PREF PIPE FLEX YOFF 


float 


Measured flexure ofl'set in z/-direction in CCD pixels 


hierarch PPAK PREF PIPE SPEC RES 


float 


Homogenized spectral resolution (FWHM) in A 


hierarch PPAK PREF PIPE SKY MEAN 


float 


Mean surface brightness (//sky) during observation 



set of galaxies to be released, a detailed characterization of the 
data quality for the existing data is needed to select and define a 
minimal useful suite of quality control (QC) parameters. 

Overall, the final data quality for CALIFA depends on a 
number of independent factors such as (i) general instrument re- 
liability and temperature stability, (ii) ambient conditions dur- 
ing observations, and (iii) the robustness of the data reduction 
pipeline. We define a set of parameters and QC flags in the areas 
of astrometric, spectroscopic and photometric characterization 
that are measured for all galaxies. We amend these parame- 
ters with four quality flags for each of the two spectral setups 
(V500/V1200; eight in total) that summarize the reduction pro- 
cess for each galaxy and indicate data quality. 

We strongly urge users to consider these QC parameters with 
regard to their science application. For example, if a science goal 
focuses primarily on kinematics, then a substantial ofl'set in ab- 
solute astrometric registration is not a limitation for the analysis 
of any one galaxy. However, if, under the same circumstance, 
the science goal requires deriving joint spatially-extended infor- 
mation from matching with the SDSS, then caution is advised. 

The values and uncertainties listed for the individual QC pa- 
rameters are intended to provide the user with a summary of the 
information content of the CALIFA data on a specific galaxy 
in order to allow for the initial assessment of whether or not a 
specific datacube can be used for a planned analysis. In the fol- 
lowing we display and comment on the distribution of the key 
flags and QC parameters while demonstrating the overall sam- 
ple properties and the nature of outliers. 



The QC parameters and QC flags can be found in Tables [5] 
and |6] for the V500 and VI 200 data, respectively, which are also 
available in electronic form on the DRl web page and accessible 
through the Virtual Observatory (see Sect. [63] for details). 



5.1. Initial reduction pipeline flags 

As an initial description to evaluate the data and reduction qual- 
ity, we provide each of the four flags per setup that come directly 
from the reduction process. Some of the extended parameters in 
the sections below will add to further validate and expand on 
these initial flags. 

Potential minor issues are noted with a value=l, while the 
standard case receives value=0. Galaxies with major issues, i.e. 
flag values =2 have been observed as part of the CALIFA survey, 
but were excluded from this data release. The reduction flags 
are: 

- FLAG_RED_0: observing conditions quality. 

Set to 1 if the night sky is brighter than jiy = 20.5 
mag arcsec"^ (V500) or jib = 20.0 mag arcsec"^ (V1200), 
respectively. There exists a correlation between limiting con- 
tinuum sensitivity of the datacubes and the night- sky bright- 
ness. For nights brighter than this limit the datacubes may 
not be as deep as targeted. 

Flag is set to 2 if observing conditions were non-photometric 
or if background of unknown origin is seen in the data. 

- FLAG_RED_R: reduction/calibration quality. 
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Set to 1 if reduction deviates from the standard procedure, 
e.g. due to a lack of directly associated continuum or arc- 
lamp exposures or a saturation in the available calibration 
data. 

Flag is set to 2 if manual intervention was required to recover 
a minimum of possible data from the raw data. 

- FLAG_RED_W: wavelength calibration quality. 

Set to 1 if the standard deviation of the pipeline-estimated 
systemic velocity in different wavelength regions is larger 
than 25 kms"^ 

Flag set to 2 if the deviation is more than 2 standard devia- 
tions, or 34kms"^ 

- FLAG_RED_V: visual inspection quality. 

Summary of a visual inspection of the data. If there is an 
obvious defect affecting a small fraction of the FoV (like the 
presence of a strong line or a imperfectly traced fiber), flag 
is set to 1 . It is set to 2 when it is evident that the data have 
low quality for any unknown reason. 

Finally this flag is set to 1, even if all previous cases have 
value=0, when a master continuum or master arc frame, in- 
stead of an individual frame, has been used during the reduc- 
tion of any of the frames used to create the final datacube. 

5.2. Astrometric accuracy and spatial resolution 

5.2.1 . Astrometric registration accuracy 

The first area of testing and evaluation of the pipeline output 
is the absolute astrometric registration of the datacube coordi- 
nate systems to the International Coordinate Reference System 
(ICRS). Astrometric registration is of central importance when 
CALIFA data are to be combined with e.g. imaging data from 
other surveys to extract spatially resolved information at certain 
locations in a galaxy. For this purpose the CALIFA pipeline uses 
a simple scheme whereby the tabulated coordinates of the galaxy 
V band photometric center, as given in Table [T] are assigned to 
the measured barycenter of the reconstructed image in the CAL- 
IFA datacubes. 

We tested independently the robustness and accuracy of this 
approach by both visual inspection of the assigned galaxy cen- 
ters in both setups, as well as in V500 matching the galaxy loca- 
tion of a reconstructed g band from CALIFA to, where available, 
corresponding SDSS images. The deviations of the CALIFA as- 
trometry from the ICRS are typically small. 

In the V500 setup the offset for most galaxies lies below 
1V4. However, three galaxies show offsets between 1V4 and 3'', 
and further two substantially larger offsets of ~ 12'' and -22''. In 
these cases, either the centers of the galaxies are not well defined 
due to dust lanes in strongly inclined system, or the galaxies have 
a bright field star near their center: 

- UGC 10650: Edge-on spiral with fuzzy center. Registration 
should be better than 3". 

- NGC 4676 A: One galaxy of The Mice with unclear center 
and a possible offset of ~1'.'5 in the astrometric registration. 

- NGC 6032: Highly inclined galaxy with prominent dust lane 
and difficult to define center. ~2" offset in the coordinates. 

- NGC 0477: Astrometric center set to nearby star; resulting 
off'set is -12". 

- NGC 3991: Center of this inclined spiral is not the brightest 
component. The correct center and hence coordinate system 
is off'set by -22". 

In the VI 200 setup we tested the astrometric co-registration 
with respect to the V500 reference frame. Generally the regis- 
tration is accurate to better than 2". Only for two galaxies we 
find a more substantial offset: 
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Fig. 12. Distribution of the seeing during the CALffA observations as 
measured by the automatic Differential Image Motion Monitor (DIMM, 
|Aceituno|2004l ). 



- NGC 4470: Galaxy with difficult to define center. Offset of 
~9" between the two setups. 

- NGC 4676A: As stated for V500 the center is not well de- 
fined for this galaxy and an offset of ~6" exists between the 
two setups. 

For this DRl we provide the pipeline registration as de- 
scribed above and do not correct offsets found through these sep- 
arate external checks. In future data releases it is planned to im- 
plement more sophisticated registration methods directly in the 
pipeline thereby reducing the number of outliers and improving 
the overall astrometric registration accuracy. We describe the as- 
trometric offsets for DRl with the flag FLAG_ASTR, given for 
the V500 and VI 200 in Tables [5] and [6| respectively. For V500 
the values 0, 1, and 2 describe offsets <1'.'4, 1V4 to 3'.'0, and 
>3'.'0. In the V1200 table values for FLAG_ASTR of and 2 re- 
fer to relative offsets between V500 and VI 200 smaller or larger 
than 3". 



5.2.2. Seeing and spatial resolution 

The achievable spatial resolution of imaging data is usually de- 
termined by the telescope aperture together with the atmospheric 
and instrumental seeing during ground-based observations. In 
the case of CALIFA, the coarse sampling of the large PPak fibers 
modified by the adopted 3 -fold dithering pattern as described in 
S12 imposes additional limitations on the final spatial resolution. 

We measure the seeing in the CALIFA data from two in- 
dependent sources. The first one, the differential image motion 
monitor (DIMM, pVceituno|2004| ), operates fully automatically at 
the Calar Alto observatory during the nighj^ The DIMM seeing 
for the DRl sample is shown in Fig. [12] and has a median value 
of 0'.'9 (FWHM). We confirm this value with measurements of 



^ Since the DIMM has tighter operational constraints (humidity lower 
than 80% and wind speed less than 12ms~^) than the 3.5m telescope, 
seeing information is not available for every CALIFA observation. 
Hence DIMM seeing values can be missing from Tables [5] and |6] but 
the overall seeing distribution is not expected to be very different. 
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Fig. 13. Left panel: V band image of 
NGC 7321, an example of a CALIFA galaxy 
with a star in the FoV to measure the PSF of the 
data. Right panel: Distribution of the FWHM 
of the CALIFA PSF as measured from 34 cubes 
with a sufficiently bright star in the field. The 
median of the distribution is 3V7; the spread 
of values is due to the flux of the underlying 
galaxy structure and undersampHng of the stel- 
lar image and hence high- value outliers are up- 
per limits in these cases. 



the width of the guide star on images taken by the Acquisition 
and Guiding (AG) Camera of the PMAS instrument ( |Roth et aL] 
[2005 ), which, as expected, has a slightly larger median value of 
1V15 (FWHM). This confirms that atmospheric and instrumental 
seeing are not the limiting factors in the spatial resolution of the 
CALIFA cubes. 

Instead, the large aperture of each fiber (2V7 diameter) in the 
PPak fiber bundle does not permit to take full advantage of the 
good seeing conditions. The final spatial resolution of the CAL- 
IFA data is set by fiber size and the dither scheme together with 
the adopted image reconstruction algorithm. The width of the 
point spread function (PSF) is directly measured from the V500 
CALIFA datacubes in 35 cases where a bright field star is in 
the FoV (see left panel of Fig. [13] for an example). The mea- 
sured FWHM distribution of the PSF (Fig. [13] right panel) has 
a well-defined median value of 3V7 (FWHM). To verify these 
measurements we simulated the expected CALIFA PSF for V 
Gaussian seeing while adopting the fiber sampling, dither pat- 
tern and image reconstruction algorithm used for the DRl. We 
obtained a FWHM for the PSF of > y as a lower limit, which 
is in agreement with our empirical measurements. 

This spatial resolution is for a large part set by the adopted 
kernel for the inverse-distance weighting scheme for the image 
reconstruction (see S12 for details). It was chosen to produce 
smooth images without obvious structure caused by the dither 
pattern. Other image com bination approaches such as drizzling 
( Fmchter & Hook||2002| ) can be used to reach closer to the in- 
trinsic spatial resolution of CALIFA data, which is significantly 
better than y due to the dithering scheme used. The choice of 
image combination procedure depends on the actual goal and 
will be an area of improvement in future CALIFA data releases. 

5.3. Wavelength calibration accuracy 

During the data reduction the spectral resolution was homoge- 
nized to reach a target FWHM of 6A (V500) and 2.3A (V1200), 
respectively, over the whole wavelength range, as explained in 
S12. To go beyond the simple flag-cut from the pipeline di- 
agnostics we tested the spectral resolution and the wavelength 
calibration, and its spectral and spatial consistency. 

The distribution of spectral resolutions for all galaxies is 
shown in Fig. 14 for both instrument setups. It was estimated by 
measuring the width of night sky lines in the error cube, which 
contains a reliable trace of the sky lines in the input data. For 
the V500 setup we find that it is centered on the target spectral 
resolution (6A). This is not the case for the VI 200 setup. For 
this setup only Hg pollution lines (Hgi/14046A and /14358A) 




2 3 4 5 6 7 

Spectral resolution (FWHM) [A] 

Fig. 14. Distribution of the final spectral resolution (FWHM) in the 
combined cubes as measured from the sky lines in the error cube. The 
vertical dashed lines indicate the spectral resolution to which the spec- 
tra were aimed to be homogenized during the data reduction. Measure- 
ments for VI 200 (blue, left) are upper limits since the night sky lines 
used are intrinsically resolved at this level. The spread is due to uncer- 
tainties in calculating the spectral resolution from only a few available 
lines. 



are available in the corresponding wavelength range. Since these 
lines are resolved at the resolution of this setup the derived em- 
pirical resolution is a conservative upper limit. 

We tested the overall wavelength calibration in two ways. 
First, we measured the centroids of the same night sky lines 
and determined their scatter across the FoV. In all cases, the 
centroids are fully consistent with zero oflfset from the nominal 
wavelength, while the scatter is consistent with pure measure- 
ment errors and the absence of any detectable systematic spatial 
variation. 

Secondly, we tested whether there exists any systematic oflf- 
sets of the wavelength calibration at diflTerent wavelengths of the 
spectra. For this purpose we modeled the recession velocity in 
the spectra of every galaxy's coadded central 10'' apertures in 
diflferent spectral regions . We did this using the pena lized Pixel- 
Fitting method (pPXF, ,Cappellari & Emsellem| 2004) to piece- 
wise fit stellar models, hence measuring central wavelengths in 
several (3-5) independent spectral bins showing strong absorp- 
tion features. The resulting velocity diflTerences between these 
bins were again always consistent with pure measurement noise 
and the absence of systematics. 
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Fig. 15. Cumulative distribution of the V band (V500, left curve) 
and B band (VI 200, right curve) night- sky surface brightness during 
the observations of the V500 and VI 200 CALIFA data, respectively. 



In summary, no spatial or spectral variation of the wave- 
length calibration is found beyond the uncertainty imposed by 
the spectral resolutions of 6 A and 2. 3 A (FWHM) for the dif- 
ferent setups - hence the random shot noise of the data fully 
dominates the accuracy of wavelength calibration. 

5.4. Spectrophotometric quality and accuracy 
5.4.1 . Sky subtraction quality 

Cumulative distributions of the V and B band night- sky sur- 
face brightness during observations are shown in Fig. [15] for 
the V500 and VI 200 CALIFA data, respectively. They were 
derived from the sky-fibers following the analysis described 
in S12. There is considerable scatter in night- sky brightness 
of about 1.5 magnitudes in both cases, indicating the highly 
variable conditions at the observatory. The typical night- sky 
brightness is yUy ~ 21.0 mag arcsec"^ for the V500, and yUg ~ 
22.1 mag arcsec"^ for the VI 200 data. While the former is 
~1 mag brighter than the mean dark night sky at the observa- 
tory 22magarcsec"^), the latter is much m ore similar to 
the mean (jufi- 22.5 mag arcsec"^), according to Sanchez et al 



(|2007a). This is expected because we have reserved the dark- 
est nights (new moon ±1 night) for the VI 200 setup, while gray 
nights are used for the more efficient V500 setup. 

We have conducted a series of tests to assess the sky subtrac- 
tion accuracy in the continuum as delivered by the data reduction 
pipeline. As the initial test we again used the Starlight code 
( Cid Fernandes et al.|2005| ) to model the central pipeline-reduced 
spectra with stellar templates and allowing for the presence of a 
dust component. Figure [16] shows the residuals of these fits as 
a fixed rest- frame wavelength stack of spectra for the 100 DRl 
galaxies. Here spectra were always integrated over a central 10'' 
aperture. 

In this view, all rest-frame features form vertical stripes 
while atmospheric and instrumental features are slanted away 
from the vertical. The main effects from imperfections in sky 
subtraction can be seen around 6000 A and 625 OA, where sky- 
line residuals at a few percent of the stellar template flux exist. 

In order to be more quantitative regarding how much the sky 
subtraction could be systematically improved, if at all, we have 



taken advantage of the ability of the pPXF package to simultane- 
ously fit sky and stellar templates to optimize the sky subtraction 
of our dat43 

We selected all sky fibers from each of the three pointings 
for a particular object and, as in the pipeline, discarded the six 
brightest fibers to ensure that we eliminate cases that are clearly 
not dominated by the sky. The remaining 30 sky fibers were used 
as sky templates for pPXF. The input galaxy spectra were se- 
lected from fibers located in a ring in the inner disk of the galaxy 
where light intrinsic to the galaxy dominates. We noticed that 
the use of fibers from the outer parts of the galaxy (i.e. more sky 
dominated) can lead to errors in the derived kinematic param- 
eters, giving the wrong sky-template plus stars combination as 
the output of the pPXF fitting. After masking emission lines to 
avoid nebular contamination we ran pPXF, which gives as output 
the best stellar combination plus the best sky combination. This 
procedure therefore provides the combination of sky fibers that 
best matches the data. For further optimization, we ran pPXF on 
a fiber-per-fiber basis to adjust this best sky solution to each of 
the fibers including the sky-dominated outer fibers. The resulting 
sky weights are rather constant over the FoV of the instrument 
with a scatter of 4%, as is illustrated in Fig.[T7] 

This method was applied to a set of objects in the DRl sam- 
ple. We then compared the sky- subtracted spectra of a point- 
ing coming from the pipeline and our method, before cube re- 
construction. We carried out this comparison by calculating the 
S/N in the continuum: for those fibers for which the method 
described above improves the sky subtraction, the average dif- 
ference in the S/N is 3.7 ± 2.8%. 

Hence, the resulting sky subtraction is, globally, not signifi- 
cantly improved compared to the pipeline sky subtraction. The 
average improvement achieved in some of the fibers is of the or- 
der of a few percent, thus nearly as the pipeline sky subtraction 
is within a few percent in S/N as accurate as the sky subtraction 
carried out using a sophisticated template fitting algorithm. 

For completeness, we note that among all lines the strongest 
residual from the sky occurs in the Hgi/14358A, Hg 1/15461 A, 
O I ^5577A, Na i D (around 5890 A), O i ^6300A and O i ^6364A 
emission lines as well as at the telluric B band absorption. 

5.4.2. Spectrophotometric accuracy 

As part of the standard CALIFA pipeline each V500 datacube 
was rescaled in the absolute flux level to match the SDSS DR7 
broad-band photometry within an aperture of 30'' diameter, 
while the V1200 data were matched to the V500 data (S12). 
A mosaiced SDSS image was created for each galaxy u sing 
SDSSmosaic, an IRAlJj package developed by S. Zibetti (I Zi-| 
|betti et al.|2009 ), which also subtracts the background from each 
SDSS scan by fitting a plane surface that allows for linear gradi- 
ent along the scan direction. All stripes of the mosaic were com- 
bined using the program SWarp ( [Bertin et al.|2002| ). Finally, we 
obtained the Galactic-extinction corrected absolute magnitudes 



^ A caveat is that while ideally the following method should be imple- 
mented at the pipeline level, at this point it is too complex and needs 
too much intervention to be implemented there, or to be carried out for 
all the sample galaxies - at the same time it does not provide optimum 
sky fitting for all the fibers. 

IRAF is distributed by the National Optical Astronomy Observato- 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Sci- 
ence Foundation. 
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rest-frame wavelength [A] 

Fig. 16. Stack of 100 DRl central-region spectra in the rest-frame, sorted by redshift. Shown are the residuals after subtraction of the best-fitting 
stellar population model with Starlight. Systematic deviations from a idealized model can be seen as vertical stripes (rest-frame mismatches, e.g. 
imperfect stellar model or emission lines), slanted stripes (observed-frame mismatches, e.g. imperfect sky model) or general blue-to-red features 
(imperfect relative flux calibration with respect to SDSS). 
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Fig. 17. Left panel: Sky weight fac- 

tor of the sky spectrum as inferred with pPXF 
across the PPak field for the second pointing of 
UGC 03944. Right panel: Distribution of sky 
weight factors for that particular pointing. 
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Fig. 18. Left panel: Distribution of the 30'' aperture photometry scale factor between the SDSS DR8 images and re-calibrated CALIFA data. We 
compare the photometry only for the g and r bands, which are both entirely covered by the V500 wavelength range. Right panel: Distribution of 
the corresponding color off'set between the SDSS DR8 images and the synthetic CALIFA broad-band images. 



with SExtractor (Bertin & Amouts'1996') using the photometric 
zero-points as provided by the SDSS DR7 photometric pipeline. 

This process to rescale fluxes for absolute re-calibration will 
not be perfect as can be seen from the slight systematic residu- 
als in the blue and red edges of the spectrum in Figure [16] To 



quantitatively cross-check the achieved accuracy in the absolute 
photometry of the CALIFA data, we use the improved photome- 



try provided by the 8th data release of SDSS fA ihara eta l.'2011) 
to independently measure any systematic oflTset. The distribu- 
tion of photometric scale factors for the g and r bands is shown 
in Fig. [18] (left panel). The mean SDSS/CALIFA g and r band 
ratios in the release are 1.06 + 0.13 and 1.00 ± 0.10, respectively, 
which shows that the absolute photometric calibration of CAL- 
IFA data is better than 15%. We attribute the mild systematic 
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offset to the improvements in the SDSS sky subtraction accu- 
racy around the low-redshift galaxies with size of >W selected 
for CALIFA. However, there are also several objects with photo- 
metric scale factors below 0.8 and above 1.2 for which we iden- 
tify bright field stars as being the most likely cause of a bias in 
the background estimates either for the SDSS or CALIFA data. 
For DRl we actually drew a line at values of 0.5 and 2.0 for 
the scale factors and removed a few objects that showed values 
beyond these limits. In the right panel of Fig. 18 we show the 
distribution in A(g - r) color diff'erence between the SDSS and 
CALIFA data. We find a systematic off'set of A(g-r)=-0.06 mag 
(median) with a scatter of only 0.05 mag. Thus, the spectropho- 
tometric accuracy across most of the covered wavelength range 
is 6% for the CALIFA data. 

For the V500 setup we add a flag for the quality of photo- 
metric calibration, FLAG_FLUXCAL. If for both SDSS g and r 
band the ratio between SDSS and CALIFA lies between 0.7 and 
1.4, we assign a flag value of 0. For at least one band behaving 
worse than this, with ratios up to 2.0 or down to 0.5, we assign 
FLAG_FLUXCAL= 1 , and we reject the data for even worse val- 
ues. 



5.5. Limiting sensitivity and signal-to-noise 

The limiting continuum sensitivity of the spectrophotometri- 
cally recalibrated SDSS spectra are presented in Fig. [19] We 
compare the 3cr continuum flux density detection limit per in- 
terpolated 1 arcsec^-spaxel for the faintest regions against the 
average S/N per spectral resolution element within in an el- 
liptical annulus of ±V' around the galaxies r band half-light 
semi-major axis (Rso)- Both values were estimated from the 
noise within the narrow wavelength ranges 44 80-45 20 A and 
5590-5680A for the V1200 and V500 setup, respectively, that 
are nearly free of stellar absorption features. As expected, 
the 3(T continuum flux density detection limij^ for the V500 
data (as low as I^o- = 1.0 x 10"^^ erg s"^ cm"^ A"^ arcsec"^ in 
the median at 5635A) is a factor of ~2 fainter than for the 



V1200 data (/s^ = 2.2 x 10"^^ erg s"^ cm"^ A"^ arcsec"^ in the 
median at 4500A) mainly because of the difl'erence in spec- 
tral resolution. These continuum sensitivities can be directly 
transformed into equivalent limiting broad band surface bright- 
nesses of 23.6 mag arcsec"^ in the r band for the V500 and 
23.4 mag arcsec"^ in the g band for the VI 200 data. The dis- 
persion in the limiting continuum sensitivity can partially be ex- 
plained by the variance in the sky brightness of the correspond- 
ing night, because both are significantly correlated, as expected. 
However, other atmospheric conditions, such as dust attenuation, 
transparency of the night, etc. are important characteristics that 
alter the achievable depth at fixed exposure times. We provide a 
flag for the limiting sensitivity, FLAG_DEPTH, which indicates 
with 0, 1,2, whether the 3cr continuum limiting sensitivity is <5, 
between 5 and 20, or >20 x 10"^^ erg s"^ cm"^ A"^ arcsec"^ . 

Another interesting quantity for practical purposes is the S/N 
of the data. Since the galaxy brightness enters as the signal, S/N 
correlates only partially with the limiting sensitivity, which is a 



^ Please note that this is a continuum flux density. For e.g. any 
emission line the spectral resolution has to be taken into account, 
i.e. 6A and 2.3A FWHM for the V500 and V1200 setups. In or- 
der to detect an unresolved emission line with a 3<t significance in 
their peak flux density requires an integrated emission line flux of 



-1x10 ^^,ergs ^cm^arcsec ^ and -0.6x10 
in the V500 and VI 200 data, respectively. 
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Fig. 19. Limiting 3cr continuum sensitivity measured at the outskirts 
of the galaxy as a function of the average continuum S/N at the half light 
radius (R50). The corresponding broad-band surface brightness limits in 
r (V500) and g (VI 200) are indicated on the right ^-axis. The limiting 
continuum sensitivity and the S/N was computed from the median sig- 
nal and noise in the wavelength region 4480-4520A and 5590-5680A 
for the V1200 and V500 data, respectively. 



measure of the noise. The mean S/N in the continuum per spaxel 
at the half-light semi-major axis (IR50) of all objects is ~14 for 
the V1200 setup, while it is -34.5 for the V500 data. Thus, we 
achieve a S/N>10 for the majority of the objects even for the 
VI 200 setup. This means that using adaptive binning scheme, a 
decent S/N can be achieved far beyond IR50. 



5.6. Overall quality judgment 

The quality of a given dataset is defined by a combination of the 
quality flags introduced above. We used a subset of these in order 
to guarantee that we deliver scientifically first-grade data in this 
data release. Minimum (maximum) threshold values have been 
defined as follows: Airmass, sec(z)<2; relative spectrophotomet- 
ric calibration relative to SDSS, 0.50 < (/califa//sdss) < 2 
for both g and r band; and a limiting sensitivity better than 



2 X 10"^^ ergs 



^ cm ^ arc sec ^ 



in the abovementioned spectral 



windows (Sect. |5.5| ). We also add the described checks on see- 
ing, wavelength calibration and spectral resolution and initial re- 
duction flags. 



Article number, page 13 of 



27 



For each object we list the most important QC parameters 
discussed in this section in Tables |5] and [6] We carried out very 
extensive vetting and description in order to deliver 100 galaxies 
with prime quality data. These listed parameters can be used to 
check the data quality of the released cubes per object in more 
detail. 



6. Access to the CALIFA DR1 data 

6.1. The CALIFA DR1 search and retrieval tool 

The public data is distributed through the CALIFA DRl web 
page (http://califa.caha.es/DRl). A simple web form in- 
terface, developed especially for this first data release, allows 
to select data of a particular target galaxy, or a subsample of ob- 
jects within some constraints on galaxy properties and observing 
conditions. Among the selection parameters we include the in- 
strument setup, galaxy coordinates, redshift, g band magnitudes, 
observing date, Hubble type, bar strength, inclination estimated 
from axis ratio, V band atmospheric attenuation, airmass, and 
relative accuracy of the SDSS/CALIFA photometric calibration. 
See W12 for a more detailed characterization of the CALIFA 
galaxies. 

If CALIFA datasets are available within the given con- 
straints, they are listed in the succeeding web page to make a 
final selection for download. The download process requests 
a target directory on the local machine to store the data, after 
the downloading option was selected. The CALIFA data are de- 
livered as fully reduced datacubes in FITS format, described in 
Sect. [3] separately for each of the two CALIFA spectral settings, 
i.e. the V500 and VI 200 setup. Each DRl dataset is uniquely 
identified by their file name, GALA^AM£^.V1 200. rscube. fits. gz 
and GALA^AM^.VSOO.rscube.fits.gz for the V1200 and V500 
setup respectively, where GALNAME is the name of the CAL- 
IFA galaxy listed in Table [T] 

6.2. CALIFA Galaxy Explorer page 

As it can be a challenge for the general user to judge a-priori 
which CALIFA dataset could be useful for their specific sci- 
entific interest, we prepared a CALIFA Galaxy Explorer page 
(hereafter CGE) that provides a presentation card for each indi- 
vidual CALIFA DRl galaxy to the scientific community. It is a 
powerful tool for a quick-look reference of the galaxy properties 
and the actual content of the CALIFA DRl data. The CGE is 
linked to each object in the output table of the CALIFA retrieval 
tool and a column with a link is added to a dedicated DRl sample 
table available on the DRl web page. 

The CGE page is structured in three diff'erent levels. In the 
first level, a false-color (gri bands) SDSS postage stamp image 
(1.5' X 1.5' in size) is presented together with a short table con- 
taining the galaxy name, galaxy equatorial coordinates (in J2000 
format) and the designated CALIFA ID. Alternative names of the 
galaxies can be found through a direct link to NED together with 
the redshift, the SDSS u, g, r, /, and z band Petrosian magnitudes, 
and the corresponding Galactic extinction in each band. 

In the second level, the user finds the CALIFA data prod- 
ucts and their content. A compact visualization of the CAL- 



IFA data (see Fig. 20 ) shows the reconstructed V band image, 
the stellar population age distribution, the gas velocity field, and 
the HQf-l-[N ii] /1/16543, 6584 emission-line 3D distribution at first 
glance. Those maps were created by a preliminar y analysis of 
the V500 datacubes with the FIT3D software tool ( Sanchez et al.| 
|2007b| ). Additionally, the B and V band images reconstructed 



from the CALIFA datacubes are shown for both the VI 200 and 
V500 setups, respectively, together with the corresponding 30'' 
aperture integrated spectra. 

In the third level, we report the result of our visual morpho- 
logical classification of the galaxy. More detailed galaxy infor- 
mation is shown in a separate table which includes additional 
parameters available through the Hyperleda database. 

6.3. Virtual Observatory services 

CALIFA data can be accessed using the Virtual Qbser- 



vatory (VO) Table Access Protocol (TAP, [Dowl er et al. 
[20TT] ) at the TAP URL http://dc.g-vo.org/tap (service 
ivo://org.gavo.dc/ system /tap/run). We provide tables de- 
scribing each cube as a single data product (califadrl. cubes) in- 
cluding the QC parameters reported in Tables [5] and |6] as well 
as tables containing all fluxes of the entire DRl by position 
and wavelength (calif adrl.fluxv500 and calif adrl.fluxv 1200) for 
the two setups separately. The table schemata can be obtained 
using the usual TAP means and tables can be easily retrieved 
and queried through the VO client TOPCAlj^ Example queries 
and links will be provided on the DRl web page to guide the 
user. The cubes can also be found in the GAVO data center's 



ivoa.obscore tables with a complete ObsCore metadata set (Tody 



et al. 2011 ) and the corresponding links are also stored in the cal- 



ifadrl. cubes table. 

Individual spectra can also be accessed using IVOA's 
Simple Spectral Access Protocol using the service 
ivo://org.gavo.dc/califa/q/s. This allows for easy access 
from within analysis programs like Splat and an example will 
be provided on the DRl web page. 

An overview of VO-accessible resources generated from 
CALIFA data is available at http://dc.g-v0.0rg/br0wse/califa/q. 

6.4. Visualization of CALIFA data 

Visualization of CALIFA data is essential for its correct treat- 
ment, analysis and eventual scientific interpretation. IFS visu- 
alization tools are very useful for this purpose, as they provide 
a 2D view at any wavelength slice or enable exploration of the 
spectra at any spaxel. There are various visualization tools avail- 
able in the community. Some are designed for a quick-look to 
the data, others allow for more interactivity, while yet other are 
developed for a particular project or data format. 

For the visualization of CALIFA data, stored in a FITS dat- 
acube format, we list here several tools as examples: 1) dsS|^ 
from version 7, can load and view 3D datacubes in multiple 
dimensions; 2) QFitsVieiAj^by Thomas Ott, is a generic FITS 
viewer program capable of handling datacubes and performing 
basic analysis operations; 3) E3e|j (Sanchez 2004), is a high- 
level IFS package which allows interactive visualization, spa- 
tial resampling, and basic analysis; 4) IFSvievI^ is a light 
and stand-alone Python script using Tklnter, pyfits and 
matplotlib, to visualize FITS datacubes. Its ultimate goal is 
to transfer all the functionality of E3D into a more flexible and 
easy to install tool; 5) CASA Viewer as part of the CASA soft- 
ware packag^^ it was initially designed for the handling and 



http : //www . star . bris . ac . uk/~mbt/topcat/| 
http : //hea- www . harvard . edu/RD/ds9 
http : //www . mpe . mpg . de/~ott/QFitsViewj 
http : //www . aip . de/Euro3D/E3D 



http : //www . caha . es/sa nchez/IFSview/ 
|http : //casa . nrao . edu 
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Fig. 20. Top panels: SDSS postage stamp images of CALIFA target galaxies, with a FoV of 90'' x 90'\ Bottom panels: Corresponding examples 
of quick-look CALIFA 3D visualization images shown on the CALIFA Galaxy Explorer page for four representative galaxies. From bottom to 
top it shows the reconstructed V band light distribution, the stellar population luminosity- weighted age distribution, the Ka gas velocity field, and 
the HQf+[Nii] AA6543, 6584 emission-line 3D distribution inferred from the CALIFA V500 datacube. The associated color bars are shown below 
from left to right. 



analysis of 3D radio data, but will be extended to display opti- 
cal/NIR datacubes in the upcoming version 4.0 (Kuntchner et al. 
in preparation). It can also read error vectors that c an be properly 
used in the analysis of the data; 6) PINGSof t[^(Rosales-Ortega 
|2011 ), is a set of IDL routines designed to visualize, manipu- 
late, and analyze IPS data regardless of the original instrument 
and spaxel size/shape. Given that it was originally developed 
for data from the PPak instrument ( Rosales-Qrtega et al.|2QlQ| ), 
it was recently updated to handle the CALIFA data format (see 
below). All those tools are also available through the CALIFA 
DRl web page. 

In the case of PINGSof t, the package is adapted to work na- 
tively with the CALIFA data, all routines automatically identify 
the different FITS extensions (HDUs) of the CALIFA format, al- 
lowing the user to spatially and spectrally explore the CALIFA 
data. In addition, PINGSof t includes routines to extract regions 
of interest by hand, within a given geometric aperture or based 
on a user-given mask, to integrate the spectra within a given re- 
gion or based on the S/N derived from continuum/emission line 
features; to convolve the data with a full set of narrow and broad- 
band filters for visualization or analysis purposes; to perform 
radial binnings with either fixed spatial bins or S/N floor; to de- 
rive cross-correlations velocity fields; to perform spatial adap- 
tive binning (e.g. Voronoi tessellation), as well as some miscel- 



http : //calif a . caha . es/pingsof t 



laneous codes for generic tasks in astronomy and spectroscopy. 
PINGSof t also includes a routine to split the diflferent CALIFA 
FITS extensions into individual FITS file£3 

In summary, there have been important efforts in the devel- 
opment of IFS visualization and manipulation tools. Given the 
number of diflTerent software options available in the literature, 
spatial and spectral visualization of the CALIFA data should not 
represent a limitation for potential users. 

7. Conclusions 

In this article we describe the data distributed with the first 
public data release of the Calar Alto Legacy Integral Field 
Area (CALIFA) survey. The CALIFA DRl provides science- 
grade and quality-checked integral-field spectroscopy data for 
100 nearby galaxies publicly distributed to the community at 
http : //calif a . caha . es/DRl The current DRl sample com- 
prises one sixth of the total foreseen sample, but it already covers 
a similar range of colors, magnitudes (see Fig.|2]), and morpho- 
logical types (see Fig. |4]). Thus, the released IFS data are rep- 
resentative of the typical spectroscopic properties of galaxies at 
each color/magnitude bin, although it lacks the proper statistical 

A similar Python routine that can split the FITS extensions into sep- 
arate files and convert the CALIFA data to CASA data format can be 
found in the CALIFA web page. 
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depth of the final sample. With this sample it will be possible 
to probe the evolutionary sequence across the color-magnitude 
diagram, even if it may not yet be possible to disentangle the 
eff'ects of diff'erent galaxy properties and their evolution at each 
bin (like active nuclei, morphology, presence of bars...). 

An improved data reduction pipeline (VI. 3c) has allowed us 
to increase the quality of the data beyond the original specifica- 
tions of the survey. The modifications with respect to the previ- 
ous pipeline (VI. 2) are described in Appendix [A[ The final FITS 
data format, including science data, propagated error vectors, 
masks, and error weighting factors was outlined in Sect|4] The 
quality of each datacube was cross-checked by a set of defined 
figures of merit, allowing us to select and distribute the data with 
the highest quality. Among the parameters cross-checked are the 
depth, wavelength calibration accuracy, spectral resolution, sky 
subtraction, and stability of the observations and reduction pro- 
cess. 

Despite its limited statistical scope, the CALIFA DRl sam- 
ple already comprises the largest (so far) IPS survey of nearby 
galaxies of all morphological types, sampling the full color- 
magnitude diagram, and covering the full optical extension of 
galaxies with a large wavelength coverage. The scientific use 
of this dataset is verified by several publications currently in 
preparation by the CALIFA collaboration on diff'erent science 
topics, including the ISM in early-type galaxies, dynamical and 
gas ionization state in merging galaxies, the stellar mass-growth 
in galaxies, dynamical state of galaxies, and oxygen abundance 
gradients. All those topics cover most of the science goals de- 
scribed in the official proposal for CALIFA survey, the CAL- 
IFA Red Boolp] but there is certainly much more to explore. 
We encourage the community to make use of the CALIFA data. 
In return we request to cite the technical CALIFA papers about 
the survey design ( [Sanchez et al.||2012a| , sample characteriza- 
tion (Walcher et al., in preparation), and this paper together with 
some recommended acknowledgement 

The amount of data distributed in DRl already comprises 
-200 000 independent spectra. While this is already a large sam- 
ple in itself, it represents only 17% of the total number of spec- 
tra foreseen at the completion of the CALIFA survey. Up to 
now we have acquired data for -200 galaxies in both setups, for 
which we will continue to perform the quality control tests im- 
plemented for this DR. It is expected that the second CALIFA 
data release will take place after 2013, when we have collected 
high-quality science-grade data of 300 galaxies in both setups. 
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Table 1. CALIFA DRl galaxies and their characteristics. 



Name 




Qf(J2000)b 


(5(J2000)b 


f 


< 




Mu - m/ 


type^ 


barf 




IC5376 


001 


00:01:19.778 


+34:31:32.40 


0.0168 


U2\ 


12.79 


3.25 


Sb 


A 


0.27 


NGC7819 


003 


00:04:24.505 


+31:28:19.22 


0.0167 


14.07 


13.61 


2.17 


Sc 


A 


0.53 


UGC00036 


007 


00:05:13.881 


+06:46:19.30 


0.0210 


14.13 


12.61 


3.41 


Sab 


AB 


0.60 


NGC0036 


010 


00:11:22.297 


+06:23:21.66 


0.0203 


13.46 


12.28 


2.83 


Sb 


B 


0.65 


UGC00312 


014 


00:31:23.921 


+08:28:00.23 


0.0145 


13.86 


13.61 


1.58 


Sd 


B 


0.35 


NGC0444 


039 


01:15:49.562 


+31:04:50.24 


0.0161 


14.55 


13.90 


2.07 


Scd 


A 


0.24 


NGC0477 


042 


01:21:20.483 


+40:29:17.33 


0.0196 


14.37 


13.20 


2.86 


Sbc 


AB 


0.66 


IC1683 


043 


01:22:38.929 


+34:26:13.65 


0.0162 


14.04 


12.75 


3.05 


Sb 


AB 


0.59 


UGC01057 


053 


01:28:53.253 


+ 13:47:37.67 


0.0212 


14.44 


13.49 


2.42 


Sc 


AB 


0.30 


NGC0776 


073 


01:59:54.525 


+23:38:39.39 


0.0164 


13.38 


12.31 


2.94 


Sb 


B 


0.69 


UGC01938 


088 


02:28:22.136 


+23:12:52.65 


0.0213 


14.50 


13.40 


2.76 


Sbc 


AB 


0.25 


NGC1056 


100 


02:42:48.312 


+28:34:26.96 


0.0052 


12.95 


11.50 


3.13 


Sa 


A 


0.57 


NGC1349 


127 


03:31:27.511 


+04:22:51.24 


0.0220 


13.39 


12.05 


3.32 


E6 


A 


0.89 


UGC03253 


146 


05:19:41.885 


+84:03:09.43 


0.0138 


13.56 


12.57 


2.77 


Sb 


B 


0.62 


NGC2410 


151 


07:35:02.261 


+32:49:19.56 


0.0156 


13.29 


12.00 


3.05 


Sb 


AB 


0.32 


UGC03995 


155 


07:44:09.127 


+29:14:50.75 


0.0159 


13.13 


11.98 


3.85 


Sb 


B 


0.46 


NGC2449 


156 


07:47:20.299 


+26:55:48.70 


0.0163 


13.48 


12.12 


3.46 


Sab 


AB 


0.50 


IC2487 


273 


09:30:09.166 


+20:05:27.04 


0.0145 


13.87 


13.02 


3.21 


Sc 


AB 


0.22 


IC0540 


274 


09:30:10.338 


+07:54:09.90 


0.0069 


14.28 


12.91 


3.16 


Sab 


AB 


0.30 


NGC2916 


277 


09:34:57.601 


+21:42:18.94 


0.0124 


13.10 


11.94 


2.80 


Sbc 


A 


0.59 


UGC05358 


306 


09:58:47.135 


+ 11:23:19.31 


0.0097 


15.21 


14.99 


1.46 


Sd 


B 


0.33 


UGC05359 


307 


09:58:51.646 


+ 19:12:53.91 


0.0283 


14.66 


13.58 


2.70 


Sb 


B 


0.37 


UGC05396 


309 


10:01:40.484 


+ 10:45:23.13 


0.0181 


14.38 


13.66 


2.71 


Sbc 


AB 


0.27 


NGC3160 


319 


10:13:55.114 


+38:50:34.53 


0.0229 


14.54 


12.98 


3.45 


Sab 


AB 


0.29 


UGC05598 


326 


10:22:14.003 


+20:35:21.87 


0.0188 


14.72 


13.68 


2.61 


Sb 


A 


0.30 


UGC05771 


341 


10:37:19.339 


+43:35:15.32 


0.0248 


14.09 


12.56 


3.40 


E6 


A 


0.71 


UGC06036 


364 


10:55:55.261 


+36:51:41.46 


0.0218 


14.07 


12.44 


3.50 


Sa 


A 


0.29 


NGC3991 


475 


11:57:30.959 


+32:20:13.28 


0.0108 


13.84 


13.37 


1.38 


Sm 


A 


0.22 


NGC4003 


479 


11:57:59.033 


+23:07:29.63 


0.0219 


14.00 


12.59 


3.27 


SOa 


B 


0.42 


UGC07012 


486 


12:02:03.145 


+29:50:52.73 


0.0102 


14.35 


14.01 


1.41 


Scd 


AB 


0.54 


NGC4185 


515 


12:13:22.192 


+28:30:39.46 


0.0130 


13.18 


12.30 


3.30 


Sbc 


AB 


0.64 


NGC4210 


518 


12:15:15.842 


+65:59:07.15 


0.0091 


13.19 


12.25 


2.84 


Sb 


B 


0.73 


IC0776 


528 


12:19:03.120 


+08:51:22.15 


0.0081 


14.77 


15.71 


0.47 


Sdm 


A 


0.56 


NGC4470 


548 


12:29:37.778 


+07:49:27.12 


0.0079 


12.81 


12.18 


1.72 


Sc 


A 


0.66 


NGC4676A 


577 


12:46:10.107 


+30:43:54.89 


0.0222 


14.71 


13.36 


2.89 


Sdm(x) 


AB 


0.28 


UGC08234 


607 


13:08:46.505 


+62:16:18.09 


0.0270 


13.54 


12.37 


2.86 


SO 


A 


0.63 


NGC5000 


608 


13:09:47.486 


+28:54:24.99 


0.0187 


13.84 


12.90 


2.70 


Sbc 


B 


0.60 


UGC08250 


609 


13:10:20.137 


+32:28:59.47 


0.0176 


15.15 


14.28 


2.26 


Sc 


A 


0.19 


UGC08267 


610 


13:11:11.334 


+43:43:34.78 


0.0242 


14.82 


13.28 


3.54 


Sb 


AB 


0.20 


UGC08733 


657 


13:48:38.994 


+43:24:44.82 


0.0078 


14.54 


14.03 


1.56 


Sdm 


B 


0.49 


IC0944 


663 


13:51:30.867 


+ 14:05:31.95 


0.0234 


13.84 


12.31 


3.73 


Sab 


A 


0.30 


UGC08778 


664 


13:52:06.668 


+38:04:01.27 


0.0108 


14.15 


12.95 


2.89 


Sb 


A 


0.21 


UGC08781 


665 


13:52:22.745 


+21:32:21.66 


0.0253 


13.82 


12.64 


3.17 


Sb 


B 


0.52 


NGC5378 


676 


13:56:51.013 


+37:47:50.05 


0.0100 


13.41 


12.36 


3.14 


Sb 


B 


0.63 


NGC5394 


680 


13:58:33.652 


+37:27:12.54 


0.0114 


14.13 


13.20 


3.12 


Sbc(x) 


B 


0.74 


NGC5406 


684 


14:00:20.119 


+38:54:55.52 


0.0180 


13.08 


11.91 


3.24 


Sb 


B 


0.88 


NGC5682 


758 


14:34:44.978 


+48:40:12.83 


0.0076 


14.38 


13.97 


1.67 


Scd 


B 


0.31 


NGC5720 


764 


14:38:33.281 


+50:48:54.87 


0.0260 


13.97 


12.85 


2.88 


Sbc 


B 


0.65 


UGC09476 


769 


14:41:32.028 


+44:30:45.97 


0.0109 


13.56 


12.85 


2.26 


Sbc 


A 


0.63 


UGC09665 


783 


15:01:32.464 


+48:19:10.92 


0.0085 


14.15 


12.89 


2.99 


Sb 


A 


0.23 


UGC09873 


797 


15:29:50.650 


+42:37:44.10 


0.0188 


15.09 


14.32 


2.60 


Sb 


A 


0.21 


UGC09892 


798 


15:32:51.947 


+41:11:29.28 


0.0189 


14.73 


13.75 


2.59 


Sbc 


A 


0.29 


NGC5966 


806 


15:35:52.108 


+39:46:08.04 


0.0151 


13.36 


11.95 


3.40 


E4 


A 


0.60 


NGC6032 


820 


16:03:01.124 


+20:57:21.32 


0.0145 


13.92 


12.85 


3.33 


Sbc 


B 


0.30 


UGC10205 


822 


16:06:40.180 


+30:05:56.65 


0.0219 


13.73 


12.36 


3.27 


SOa 


A 


0.58 


NGC6063 


823 


16:07:12.993 


+07:58:44.36 


0.0095 


13.57 


12.85 


2.54 


Sbc 


A 


0.60 


IC1199 


824 


16:10:34.346 


+ 10:02:25.32 


0.0158 


13.80 


12.60 


2.79 


Sb 


AB 


0.41 


NGC6081 


826 


16:12:56.857 


+09:52:01.57 


0.0171 


13.64 


12.10 


3.55 


SOa 


A 


0.46 


UGC10331 


828 


16:17:21.123 


+59:19:12.46 


0.0152 


14.38 


13.85 


1.84 


Sc(x) 


AB 


0.26 


NGC6125 


829 


16:19:11.535 


+57:59:02.89 


0.0154 


12.94 


11.50 


3.38 


El 


A 


0.91 
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Table 1. continued. 



Name 


ID'' 


Qf(J2000)^ 


(5(J2000)^ 








Mu - 


type^ 


bar 




NGC6146 


832 


16:25:10.330 


+40:53:34.32 


0.0292 


13.32 


11.87 


3.38 


E5 


A 


0.77 


NGC6154 


833 


16:25:30.483 


+49:50:24.93 


0.0199 


13.73 


12.77 


2.98 


Sab 


B 


0.65 


NGC6150 


835 


16:25:49.965 


+40:29:19.41 


0.0292 


14.00 


12.50 


3.52 


E7 


A 


0.48 


UGC10384 


837 


16:26:46.684 


+ 11:34:48.96 


0.0165 


14.54 


13.31 


2.42 


Sb 


A 


0.22 


NGC6173 


840 


16:29:44.875 


+40:48:41.96 


0.0294 


13.20 


11.78 


3.48 


E6 


A 


0.65 


UGC 10650 


843 


17:00:14.582 


+23:06:22.83 


0.0099 


15.22 


22.71 


-5.68 


Scd 


A 


0.20 


UGC 10693 


845 


17:04:53.020 


+41:51:55.76 


0.0280 


13.50 


12.15 


3.47 


E7 


AB 


0.68 


UGC 10695 


846 


17:05:05.573 


+43:02:35.36 


0.0280 


14.06 


12.67 


3.29 


E5 


A 


0.67 


UGC10710 


847 


17:06:52.521 


+43:07:19.96 


0.0280 


14.38 


13.02 


3.30 


Sb 


A 


0.25 


NGC6310 


848 


17:07:57.480 


+60:59:24.56 


0.0114 


13.52 


12.23 


3.29 


Sb 


A 


0.22 


NGC6314 


850 


17:12:38.715 


+23:16:12.29 


0.0221 


13.55 


12.23 


3.09 


Sab 


A 


0.51 


NGC6338 


851 


17:15:22.976 


+57:24:40.28 


0.0274 


13.42 


11.90 


3.82 


E5 


A 


0.66 


UGC 10796 


852 


17:16:47.724 


+61:55:12.43 


0.0102 


14.57 


14.53 


1.43 


Scd 


AB 


0.49 


UGC10811 


854 


17:18:43.725 


+58:08:06.43 


0.0291 


14.49 


13.20 


2.98 


Sb 


B 


0.42 


IC1256 


856 


17:23:47.285 


+26:29:11.48 


0.0159 


13.87 


12.89 


2.57 


Sb 


AB 


0.59 


NGC6394 


857 


17:30:21.423 


+59:38:23.61 


0.0284 


14.48 


13.19 


3.08 


Sbc 


B 


0.29 


UGC 10905 


858 


17:34:06.437 


+25:20:38.28 


0.0265 


13.71 


12.23 


3.42 


SOa 


A 


0.53 


NGC6411 


859 


17:35:32.849 


+60:48:48.25 


0.0123 


12.78 


11.54 


3.53 


E4 


A 


0.68 


NGC6427 


860 


17:43:38.598 


+25:29:38.17 


0.0108 


13.34 


11.90 


3.18 


SO 


AB 


0.61 


NGC6497 


863 


17:51:17.966 


+59:28:15.14 


0.0105 


13.59 


12.34 


3.57 


Sab 


B 


0.65 


NGC6515 


864 


17:57:25.195 


+50:43:41.24 


0.0228 


13.54 


12.22 


3.29 


E3 


A 


0.78 


UGC 11228 


865 


18:24:46.259 


+41:29:33.85 


0.0194 


13.81 


12.36 


3.40 


SO 


B 


0.57 


UGC 11262 


866 


18:30:35.698 


+42:41:33.70 


0.0186 


15.01 


14.59 


2.19 


Sc 


A 


0.39 


NGC6762 


867 


19:05:37.089 


+63:56:02.79 


0.0098 


13.87 


12.49 


3.13 


Sab 


A 


0.49 


UGC 11 649 


872 


20:55:27.619 


-01:13:30.87 


0.0127 


13.57 


12.44 


3.26 


Sab 


B 


0.88 


NGC7025 


874 


21:07:47.336 


+ 16:20:09.22 


0.0166 


12.87 


11.38 


3.67 


SOa 


A 


0.72 


UGC11717 


877 


21:18:35.412 


+ 19:43:07.39 


0.0212 


14.36 


12.81 


3.72 


Sab 


A 


0.47 


MCG-01-54-016 


878 


21:25:59.970 


-03:48:32.26 


0.0098 


14.92 


14.67 


1.52 


Scd 


A 


0.13 


NGC7194 


881 


22:03:30.937 


+ 12:38:12.41 


0.0272 


13.66 


12.15 


3.45 


E3 


A 


0.79 


UGC11958 


883 


22:14:46.882 


+ 13:50:27.13 


0.0262 


14.15 


12.84 


4.59 


SO(x) 


A 


0.74 


NGC7321 


887 


22:36:28.022 


+21:37:18.35 


0.0238 


13.45 


12.35 


2.94 


Sbc 


B 


0.69 


UGC12127 


888 


22:38:29.421 


+35:19:46.89 


0.0275 


13.59 


12.19 


3.86 


El 


A 


0.85 


UGC12185 


890 


22:47:25.063 


+31:22:24.67 


0.0222 


14.12 


12.99 


3.07 


Sb 


B 


0.47 


NGC74JdB 


893 


22:57:57.546 


+26:09:00.01 


0.0246 


13.40 


1 1.84 


3.92 


E2(x) 


A 


0.90 


NGC7466 


896 


23:02:03.464 


+27:03:09.34 


0.0251 


14.09 


12.92 


2.93 


Sbc 


A 


0.53 


NGC7549 


901 


23:15:17.270 


+ 19:02:30.43 


0.0157 


13.59 


12.53 


2.88 


Sbc 


B 


0.75 


NGC7563 


902 


23:15:55.927 


+ 13:11:46.03 


0.0143 


13.30 


11.83 


3.62 


Sa 


B 


0.68 


NGC7591 


904 


23:18:16.259 


+06:35:08.86 


0.0165 


13.44 


12.18 


2.90 


Sbc 


B 


0.59 


UGC12864 


935 


23:57:23.920 


+30:59:31.45 


0.0156 


14.35 


13.95 


2.51 


Sc 


B 


0.32 


NGC4676Bh 


939 


12:46:11.235 


+30:43:21.87 


0.0218 


15.03 


13.68 


3.94 


...(X) 




0.82 



Notes. CALIFA unique ID number for the galaxy. Equatorial coordinates of the galaxies as provided by NED. Redshift of the galaxies 
based on SDSS DR7 spectra or complemented with SIMBAD information if SDSS spectra are not available. Petrosian magnitudes as given 
by SDSS DR7 database corrected for Galactic extinction. Morphological type from our own visual classification (see W12 for details). Bar 
strength of the galaxy as an additional outcome of our visual classification. Ratio between the semi-minor and semi-major axis based on a 
detailed re-analysis of the SDSS images (see W12 for details). A visual morphological classification of this particular galaxy NGC 467 6B is 
missing. 
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Table 5. CALIFA DRl quality control parameters for the V500 data 



ID 


airmass^ 




Ay' 


seeing^ 


flags(R)^ 




1 ^CALIFA \g 
\ ^SDSS / 


j ''CALIFA \h 
\ ''SDSS / 




hJ 


flags(QC)^ 


001 


1.02 + 0.01 


21.4 


0.18 




0000 


5.95 + 0.33 


1.01 


0.94 


35+6 


0.8+0.2 


0000 


003 


1.01+0.01 


20.8 


0.23 


0.8 + 0.1 


0000 


6.46 + 0.98 


1.01 


1.05 


18+5 


1.0+0.3 


0000 


007 


1.34 + 0.05 


20.5 


0.17 




0000 


6.11+0.41 


0.92 


0.95 


32+5 


1.1+0.3 


0000 


010 


1.31+0.05 


21.0 


0.22 




0000 


5.70 + 0.19 


1.05 


0.98 


27+7 


1.0+0.2 


0000 


014 


1.26 + 0.04 


21.0 


0.21 




0000 


5.72 + 0.33 


1.43 


1.31 


28+7 


0.6+0.2 


0010 


039 


1.03 + 0.01 


20.8 


0.23 


0.9 + 0.1 


0010 


6.02 + 0.18 


1.04 


1.05 


19+4 


0.9+0.3 


0000 


042 


1.10 + 0.03 


21.0 


0.13 




0000 


6.77 + 0.57 


1.06 


1.00 


17+3 


1.1+0.3 


0200 


043 


1.01+0.01 


21.3 


0.19 




0010 


6.12 + 0.22 


0.99 


0.91 


37+6 


0.8+0.2 


0000 


053 


1.19 + 0.04 


20.7 


0.17 


0.8 + 0.1 


0000 


6.35 + 0.29 


1.04 


1.05 


32+7 


1.0+0.2 


0000 


073 


1.04 + 0.01 


21.1 


0.20 




0000 


6.47 + 0.45 


1.03 


0.95 


24+4 


1.5+0.3 


0000 


088 


1.13 + 0.03 


20.8 


0.17 


1.0 + 0.4 


0000 


5.65 + 0.87 


1.04 


0.99 


27+5 


1.2+0.3 


0000 


100 


1.37 + 0.26 


20.7 


0.14 




0000 


6.39 + 0.37 


1.00 


0.99 


36+7 


1.9+0.5 


0000 


127 


1.22 + 0.02 


20.6 


0.28 


0.9 + 0.1 


0000 


5.90 + 0.32 


1.00 


1.00 


19+3 


2.0+0.5 


0000 


146 


1.47 + 0.00 


20.7 


0.12 




0000 


4.86 + 0.79 


0.97 


0.94 


21+4 


1.7+0.4 


0000 


151 


1.07 + 0.02 


21.0 


0.15 




0000 


7.01+0.16 


1.07 


1.03 


26+7 


1.6+0.4 


0000 


155 


1.03 + 0.01 


21.1 


0.16 




0100 


6.34 + 0.28 


1.05 


1.02 


18+5 


1.6+0.5 


0000 


156 


1.12 + 0.04 


21.0 


0.16 




0100 


6.33 + 0.36 


0.99 


0.94 


22+14 


1.4+0.4 


0000 


273 


1.06 + 0.01 


21.1 


0.14 




0000 


6.67 + 0.19 


1.05 


1.00 


27+5 


1.0+0.3 


0000 


274 


1.18 + 0.02 


20.9 


0.17 




0000 


5.64 + 0.29 


1.03 


0.95 


34+7 


0.9+0.2 


0000 


277 


1.60 + 0.12 


20.3 


0.33 




1000 


6.69 + 0.47 


1.07 


1.12 


21+4 


1.8+0.6 


1000 


306 


1.12 + 0.01 


20.9 


0.16 


1.0 + 0.1 


0000 


5.54 + 0.06 


1.05 


0.99 


13+5 


0.8+0.3 


0000 


307 


1.06 + 0.01 


21.2 


0.16 




0100 


3.84 + 0.71 


1.07 


1.04 


18+4 


1.2+0.3 


0000 


309 


1.27 + 0.05 


21.0 


0.15 




0000 


7.44 + 0.30 


1.11 


1.04 


15+4 


1.6+0.4 


0000 


319 


1.11+0.03 


20.5 


0.56 




0000 


5.82 + 0.39 


1.01 


1.07 


25+6 


1.4+0.4 


0000 


326 


1.05 + 0.01 


20.8 


0.17 




0000 


5.72 + 0.09 


1.01 


0.96 


31+7 


0.8+0.2 


0000 


341 


0.78 + 0.68 


14.0 


0.11 


1.5 + 0.2 


1000 


6.24 + 0.41 


0.99 


0.97 


30+6 


1.0+0.3 


0000 


364 


1.33 + 0.07 


20.8 


0.12 




0000 


6.30 + 0.26 


1.00 


0.95 


39+9 


1.0+0.2 


0000 


475 


1.35 + 0.07 


20.8 


0.12 




0000 


6.42 + 0.56 


1.82 


1.66 


29+17 


1.0+0.2 


0210 


479 


1.14 + 0.03 


20.8 


0.17 




0000 


5.62 + 0.12 


1.01 


0.96 


32+6 


0.8+0.2 


0000 


486 


0.70 + 0.61 


14.0 


0.10 


1.2 + 0.7 


0000 


5.88 + 0.48 


1.03 


0.98 


27+7 


0.7+0.2 


0000 


515 


1.01+0.00 


20.9 


0.14 




0000 


6.22 + 0.24 


1.12 


1.06 


21+4 


1.7+0.3 


0000 


518 


1.16 + 0.01 


20.9 


0.15 




0000 


6.52 + 0.30 


1.03 


1.00 


23+4 


1.6+0.4 


0000 


528 


1.14 + 0.01 


20.8 


0.17 




0001 


6.30 + 0.38 


1.17 


1.16 


9+3 


0.9+0.3 


0000 


548 


1.16 + 0.01 


20.8 


0.19 


1.2 + 0.1 


0000 


6.04 + 0.50 


1.02 


0.98 


53+10 


0.8+0.4 


0000 


577 


1.02 + 0.01 


21.0 


0.15 




0000 


6.12 + 0.41 


1.07 


1.00 


17+11 


0.8+0.2 


0100 


607 


1.22 + 0.03 


21.3 


0.16 




0100 


6.26 + 0.17 


1.00 


0.94 


50+11 


1.1+0.3 


0000 


608 


1.31 + 0.06 


21.1 


0.17 




0000 


7.20 + 0.24 


1.12 


1.08 


17+4 


1.3+0.3 


0000 


609 


1.07 + 0.02 


21.3 


0.14 




0100 


4.56 + 1.09 


1.07 


1.00 


21+6 


1.0+0.2 


0000 


610 


1.11+0.03 


21.2 


0.17 




0000 


6.99 + 0.23 


1.12 


1.06 


26+7 


1.2+0.3 


0000 


657 


1.02 + 0.01 


21.2 


0.19 


1.1+0.1 


0000 


6.44 + 0.84 


1.08 


1.01 


15+4 


0.8+0.3 


0000 


663 


1.20 + 0.04 


20.9 


0.17 




0000 


5.61 + 0.42 


1.06 


0.99 


29+7 


1.0+0.3 


0000 


664 


1.04 + 0.02 


21.2 


0.17 


1.1+0.1 


0000 


6.46 + 0.96 


1.02 


0.97 


40+6 


0.8+0.2 


0000 


665 


1.06 + 0.01 


21.1 


0.17 




0000 


6.29 + 0.44 


1.03 


0.96 


22+3 


1.1+0.4 


0000 


676 


1.03 + 0.03 


21.0 


0.15 




1100 


6.62 + 1.04 


1.02 


0.97 


17+4 


1.5+0.3 


0000 


680 


1.07 + 0.02 


21.0 


0.17 




0000 


5.57 + 0.13 


1.00 


0.94 


19+9 


0.7+0.2 


0000 


684 


1.05 + 0.02 


21.0 


0.17 




0000 


5.48 + 0.11 


1.03 


0.98 


33+6 


1.5+0.3 


0000 


758 


1.02 + 0.00 


21.1 


0.24 


1.1+0.1 


0000 


7.11 + 1.14 


1.13 


1.08 


23+6 


0.7+0.3 


0000 


764 


1.29 + 0.05 


21.1 


0.18 


1.0 + 0.1 


0000 


5.43 + 0.31 


1.04 


0.94 


24+4 


0.9+0.2 


0000 


769 


1.02 + 0.01 


21.2 


0.17 




0000 


7.20 + 0.62 


1.05 


1.00 


22+4 


1.2+0.3 


0000 


783 


1.08 + 0.02 


21.3 


0.16 




0000 


5.76 + 0.91 


1.08 


1.07 


35+8 


0.9+0.2 


0000 


797 


1.05 + 0.07 


21.0 


0.15 




0000 


6.79 + 0.14 


1.01 


0.98 


23+5 


0.7+0.2 


0000 


798 


1.05 + 0.02 


21.3 


0.21 




0000 


5.68 + 0.46 


1.03 


1.00 


30+6 


0.6+0.2 


0000 


806 


1.00 + 0.00 


21.1 


0.25 


1.1+0.1 


0000 


6.33 + 0.82 


1.15 


1.08 


22+5 


1.0+0.2 


0000 


820 


1.07 + 0.02 


21.2 


0.17 




0100 


4.49 + 1.45 


1.02 


1.00 


11+3 


2.1+0.5 


0100 


822 


1.03 + 0.01 


21.2 


0.17 




0000 


7.53 + 0.17 


1.32 


1.22 


19+5 


1.4+0.3 


0000 


823 


1.15 + 0.00 


20.9 


0.19 




0000 


5.40 + 0.08 


1.09 


0.99 


26+5 


0.9+0.3 


0000 


824 


1.14 + 0.01 


20.9 


0.17 




0010 


5.60 + 0.12 


1.42 


1.13 


38+7 


0.8+0.3 


0010 


826 


1.18 + 0.03 


21.1 


0.17 




0000 


6.54 + 0.38 


1.01 


0.94 


33+7 


1.1+0.3 


0000 


828 


1.21+0.03 


21.3 


0.26 


0.9 + 0.1 


0000 


5.70 + 0.07 


1.43 


1.28 


33+9 


0.6+0.2 


0010 


829 


1.09 + 0.01 


21.2 


0.18 




0000 


5.85 + 0.21 


1.12 


1.05 


41+7 


1.3+0.2 


0000 
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Table 5. continued. 



ID 


airmass^ 


b 


A C 
Ay 


seeing^ 


nags(K) 


j: f 
Oa 


/ 9'CALIFA \g 
V ^SDSS / 


/ '"CALIFA \h 
\ '"SDSS / 


S/N(i?5o) 




ilags(QC) 


832 


1.01 +0.01 


21.3 


0.17 




0000 


7.29 + 0.68 


1.01 


0.99 


44+10 


1.2+0.4 


0000 


833 


1.12 + 0.03 


21.2 


0.18 


0.9 + 0.1 


0000 


6.16 + 0.15 


1.05 


0.94 


22+4 


1.1+0.2 


0000 


835 


1.00 + 0.00 


21.1 


0.17 


1.1 +0.2 


0000 


6.09 + 0.50 


1.01 


0.97 


43+8 


0.9+0.2 


0000 


837 


1.11 +0.00 


21.1 


0.17 


0.9 + 0.1 


0000 


5.82 + 0.12 


1.13 


1.08 


49+7 


0.7+0.2 


0000 


840 


1.18 + 0.04 


21.0 


0.28 


0.9 + 0.1 


0000 


6.26 + 0.44 


1.08 


0.99 


32+6 


1.3+0.3 


0000 


843 


1.45 + 0.09 


21.0 


0.23 


0.8 + 0.1 


0000 


6.04 + 0.92 


1.12 


0.97 


18+7 


1.0+0.3 


0100 


845 


1.03 + 0.02 


21.4 


0.19 




0000 


5.64 + 0.13 


1.08 


1.02 


40+9 


0.8+0.3 


0000 


846 


1.06 + 0.02 


21.0 


0.39 


0.8 + 0.1 


0000 


5.72 + 0.11 


1.07 


1.00 


22+4 


1.1+0.3 


0000 


847 


1.31 +0.09 


21.3 


0.15 


0.7 + 0.1 


0100 


5.41 +0.56 


0.96 


0.81 


28+11 


0.9+0.3 


0000 


848 


1.10 + 0.00 


21.0 


0.38 


0.9 + 0.1 


0000 


5.89 + 0.21 


1.06 


1.00 


36+5 


0.9+0.3 


0000 


850 


1.05 + 0.01 


21.0 


0.29 




0000 


6.35 + 0.62 


1.02 


0.99 


37+8 


1.0+0.4 


0000 


851 


1.32 + 0.04 


20.5 


0.41 


0.7 + 0.1 


0000 


5.80 + 0.08 


1.08 


1.05 


17+4 


1.6+0.4 


0000 


852 


1.12 + 0.01 


21.0 


0.36 


0.9 + 0.1 


0000 


5.77 + 0.07 


1.04 


1.01 


12+5 


0.9+0.3 


0000 


854 


1.18 + 0.03 


20.7 


0.42 


0.9 + 0.1 


0000 


5.32 + 0.90 


1.05 


0.98 


24+6 


1.1+0.3 


0000 


856 


1.05 + 0.02 


21.2 


0.17 




0000 


5.53 + 0.14 


1.06 


0.98 


28+5 


0.9+0.2 


0000 


857 


1.08 + 0.00 


21.5 


0.16 


0.9 + 0.1 


0000 


6.00 + 0.13 


1.14 


1.03 


40+7 


0.7+0.2 


0000 


858 


1.10 + 0.03 


21.2 


0.15 


1.3 + 0.2 


0000 


6.07 + 0.32 


1.06 


0.97 


36+6 


1.1+0.4 


0000 


859 


1.16 + 0.02 


21.2 


0.15 


0.7 + 0.1 


0000 


5.98 + 0.47 


1.04 


0.99 


46+7 


1.1+0.3 


0000 


860 


1.04 + 0.02 


21.2 


0.29 


1.0 + 0.1 


0000 


6.21 + 0.79 


1.00 


0.84 


44+5 


1.0+0.2 


0000 


863 


1.09 + 0.01 


21.4 


0.16 


0.7 + 0.1 


0000 


5.78 + 0.15 


0.99 


0.94 


33+6 


0.8+0.2 


0000 


864 


1.06 + 0.01 


20.9 


0.28 




0100 


5.08 + 0.44 


1.05 


0.96 


25+6 


1.4+0.3 


0000 


865 


1.01 +0.01 


21.2 


0.25 




0100 


5.75 + 0.44 


0.96 


0.90 


38+6 


0.8+0.2 


0000 


866 


1.04 + 0.02 


21.3 


0.21 


1.1 +0.1 


0111 


5.36 + 0.40 


0.98 


0.88 


16+4 


0.7+0.2 


0000 


867 


1.13 + 0.01 


21.2 


0.21 


1.1 +0.1 


0000 


5.69 + 0.20 


0.92 


0.90 


43+5 


0.8+0.3 


0000 


872 


1.28 + 0.01 


20.9 


0.15 


0.6 + 0.1 


0000 


4.93 + 0.75 


1.03 


1.00 


22+4 


1.1+0.5 


0000 


874 


1.10 + 0.02 


21.0 


0.22 




0000 


5.74 + 0.51 


0.93 


0.88 


33+6 


2.4+0.1 


0000 


877 


1.24 + 0.05 


21.0 


0.22 




0100 


5.13 + 0.41 


0.77 


0.80 


26+7 


1.1+0.3 


0000 


878 


1.33 + 0.01 


20.9 


0.22 


0.9 + 0.1 


0000 


4.00 + 0.90 


1.00 


0.95 


10+12 


0.9+0.3 


0000 


881 


1.10 + 0.00 


21.1 


0.15 


0.8 + 0.1 


0000 


5.58 + 0.32 


1.11 


1.03 


45+8 


1.0+0.2 


0000 


883 


1.10 + 0.01 


21.2 


0.17 


1.0 + 0.2 


0000 


6.11 +0.21 


0.78 


0.70 


12+3 


1.6+0.5 


0000 


887 


1.08 + 0.02 


21.3 


0.17 


0.9 + 0.1 


0000 


6.43 + 0.59 


1.06 


0.98 


38+6 


0.8+0.3 


0000 


888 


1.01 +0.01 


20.9 


0.22 




0000 


5.25 + 0.41 


1.10 


1.03 


29+7 


1.3+0.3 


0000 


890 


1.08 + 0.02 


21.0 


0.22 




0000 


4.76 + 0.71 


1.04 


1.00 


35+6 


0.9+0.3 


0000 


893 


1.35 + 0.07 


20.5 


0.20 


1.0 + 0.1 


1100 


5.96 + 0.28 


1.18 


1.04 


20+8 


2.0+0.5 


0000 


896 


1.36 + 0.08 


20.8 


0.18 


0.9 + 0.1 


0000 


6.02 + 0.14 


1.00 


0.96 


27+6 


1.0+0.3 


0000 


901 


1.10 + 0.03 


21.3 


0.14 


1.0 + 0.1 


0100 


5.54 + 0.26 


1.02 


0.95 


35+9 


1.1+0.3 


0000 


902 


1.10 + 0.00 


21.2 


0.17 


0.8 + 0.1 


0000 


5.96 + 0.13 


1.01 


0.93 


41+6 


1.0+0.3 


0000 


904 


1.20 + 0.04 


21.0 


0.17 


0.9 + 0.1 


0000 


6.10 + 0.10 


1.06 


0.99 


34+9 


1.1+0.3 


0000 


935 


1.04 + 0.02 


21.0 


0.17 


1.1+0.1 


0000 


5.79 + 0.32 


1.02 


1.02 


16+6 


0.8+0.3 


0000 


939 


1.11+0.04 


20.8 


0.15 




0000 


5.96 + 0.22 


1.00 


0.92 


17+6 


1.3+0.4 


0000 



Notes. We describe the meaning of each column including the identifier of each column in the electronic table available on the DRl web 
page. Average airmass (OBS_AIR_MEAN) and rms (OBS_AIR_RMS) of the observations for the frames used to create the consid- 
ered datacube. Average night-sky surface brightness (OBS_SKY_MAG) in the V band during the observations in units of magarcsec"^. 

Average night-sky attenuation (OBS_EXT_MEAN) in the V band during the observations in magnitudes. Average natural seeing 
(OBS_SEE_RD_AVE) in the V-band during the observations in arcsec (FWHM). Reduction qruality flags, combining the four individual column 
flags (FLAG_RED_0, FLAG_RED_R, FLAG_RED_W, FLAG_RED_V) as described in Sect. |5] Average spectral resolution (CAL_SP_RES) 
and error (CAL_SP_RES_ERR) in A (FWHM), measured by fitting the night-sky emission lines with single Gaussian functions. Ratio between 
the SDSS g band flux derived from the datacube and the one derived from the SDSS images for a 30'' -diameter aperture (CAL_DR8_RATI0_G). 

Ratio between the SDSS r band flux derived from the datacube and the one derived from the SDSS images for a 30'' -diameter aperture 
(CAL_DR8_RATI0_R). Average signal-to-noise ratio (CAL_SN_MEAN_WIN) and error (CAL_SN_ERR_WIN) estimated for the full wave- 
length range at one half light radius from the center. Average flux (CAL_LS_WIN) and error (CAL_LS_ERR_WIN) at the 3cr continuum 
detection limit in units of 10~^^ erg s~^ cm~^A~^ Quality control flags, combining the four individual column flags (FLAG_AIR, FLAG_ASTR, 
FLAG_FLUXCAL, FLAG_DEPTH) as described in Sect. [5] 
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Table 6. CALIFA DRl quality control parameters for the VI 200 data 



ID 


airmass 


}^B,sky 


Av 


seeing 


flags(R) 




S/N(i?5o) 


ha 


flags(QC) 


001 


1.01+0.02 


22.4 


0.22 




0000 


3.05 + 0.44 


13+4 


2.1+0.8 


000 


003 


1.07 + 0.06 


22.0 


0.17 




0010 


2.77 + 0.17 


9+3 


2.0+0.7 


000 


007 


1.18 + 0.03 


21.7 


0.14 




0000 


2.52 + 0.13 


12+3 


3.4+1.9 


000 


010 


1.27 + 0.09 


22.2 


0.17 




0000 


2.56 + 0.15 


8+3 


3.2+0.7 


000 


014 


1.25 + 0.08 


22.2 


0.19 




0010 


2.88 + 0.28 


11+3 


1.5+0.6 


000 


039 


1.08 + 0.05 


22.1 


0.21 


1.1 + 0.2 


0000 


3.10 + 0.14 


12+3 


1.8+0.4 


000 


042 


1.21+0.10 


22.3 


0.15 


0.9 + 0.1 


0100 


3.91 + 1.11 


11+3 


1.8+0.5 


000 


043 


1.01+0.01 


22.3 


0.27 




0000 


2.89 + 0.34 


14+4 


1.8+0.6 


000 


053 


1.36 + 0.20 


21.9 


0.16 




0100 


2.62 + 0.32 


11+2 


3.0+0.8 


000 


073 


1.07 + 0.04 


22.2 


0.26 




0110 


2.98 + 0.24 


8+2 


3.4+1.1 


000 


088 


1.34 + 0.15 


21.4 


0.21 




1000 


3.51+0.44 


10+3 


4.4+2.1 


000 


100 


1.40 + 0.16 


21.8 


0.15 


1.1 + 0.1 


0100 


3.09 + 0.69 


20+5 


2.8+0.6 


000 


127 


1.39 + 0.34 


21.7 


0.15 




1000 


3.04 + 0.14 


7+2 


3.9+1.5 


000 


146 


1.48 + 0.01 


22.2 


0.14 


1.3 + 0.2 


0100 


3.04 + 0.13 


13+3 


2.1+0.6 


000 


151 


1.57 + 0.22 


21.4 


0.17 




0000 


2.68 + 0.25 


13+4 


2.8+0.8 


000 


155 


1.08 + 0.05 


22.1 


0.13 


0.9 + 0.1 


0000 


2.97 + 0.23 


6+1 


4.1+1.0 


000 


156 


1.07 + 0.05 


22.3 


0.14 




0000 


3.68 + 0.75 


9+6 


3.5+0.6 


000 


273 


1.06 + 0.02 


22.2 


0.14 


1.1 + 0.1 


0100 


2.96 + 0.24 


15+3 


1.6+0.5 


000 


274 


1.33 + 0.11 


22.0 


0.17 




0000 


2.66 + 0.16 


12+4 


3.1+0.8 


000 


277 


1.08 + 0.04 


21.7 


0.15 




0100 


3.05 + 0.37 


15+2 


3.4+1.0 


000 


306 


1.12 + 0.01 


22.1 


0.17 


1.6 + 0.4 


0010 


2.93 + 0.17 


7+2 


1.8+0.5 


000 


307 


1.18 + 0.08 


22.2 


0.15 


0.8 + 0.1 


0000 


2.44 + 0.18 


6+2 


3.8+1.3 


000 


309 


1.23 + 0.08 


22.1 


0.14 




0000 


3.25 + 0.20 


8+2 


3.1+0.5 


000 


319 


1.07 + 0.05 


22.2 


0.15 


1.1+0.1 


0000 


2.41 + 0.09 


14+4 


1.7+0.6 


000 


326 


1.29 + 0.06 


21.9 


0.17 




0100 


2.69 + 0.22 


7+3 


3.9+1.5 


000 


341 


1.23 + 0.10 


22.0 


0.16 


1.2 + 0.1 


0000 


2.76 + 0.19 


13+4 


1.8+0.6 


000 


364 


1.02 + 0.02 


22.1 


0.17 


0.9 + 0.1 


0000 


2.84 + 0.18 


18+4 


1.8+0.6 


000 


475 


1.11+0.07 


22.1 


0.17 




0000 


2.95 + 0.30 


20+14 


2.0+0.6 


000 


479 


1.16 + 0.04 


22.3 


0.17 




0110 


2.58 + 0.10 


8+3 


3.9+1.3 


000 


486 


1.19 + 0.09 


22.1 


0.15 


0.8 + 0.1 


0000 


3.13 + 0.22 


10+5 


2.2+0.6 


000 


515 


1.18 + 0.09 


22.3 


0.15 


0.9 + 0.1 


0000 


2.82 + 0.29 


6+1 


3.7+0.9 


000 


518 


1.18 + 0.03 


22.6 


0.13 




0000 


3.22 + 0.11 


15+4 


2.1+0.6 


000 


528 


1.48 + 0.08 


22.0 


0.17 




0010 


2.37 + 0.02 


3+1 


3.4+1.1 


000 


548 


1.19 + 0.05 


21.9 


0.18 


0.8 + 0.1 


0000 


3.08 + 0.31 


19+7 


2.9+0.7 


020 


577 


1.06 + 0.04 


22.2 


0.22 




0000 


3.05 + 0.13 


8+7 


1.5+0.4 


020 


607 


1.24 + 0.07 


22.3 


0.14 




0000 


2.89 + 0.34 


23+4 


2.5+0.6 


000 


608 


1.22 + 0.11 


21.9 


0.16 




0000 


2.60 + 0.27 


9+3 


2.1+0.8 


000 


609 


1.07 + 0.05 


22.4 


0.15 


0.8 + 0.1 


0000 


2.75 + 0.20 


8+2 


2.8+0.5 


000 


610 


1.19 + 0.08 


22.5 


0.17 




0000 


2.38 + 0.05 


10+4 


2.6+0.8 


000 


657 


1.07 + 0.05 


22.4 


0.19 


1.5 + 0.2 


0000 


2.84 + 0.11 


8+2 


1.6+0.4 


000 


663 


1.10 + 0.01 


22.0 


0.17 




0000 


3.56 + 0.27 


12+3 


2.2+0.6 


000 


664 


1.08 + 0.05 


22.0 


0.26 


1.2 + 0.1 


0000 


2.82 + 0.21 


19+3 


2.0+0.6 


000 


665 


1.18 + 0.11 


22.1 


0.10 




0000 


2.69 + 0.13 


6+3 


4.0+1.1 


000 


676 


1.01+0.01 


22.5 


0.17 


0.8 + 0.1 


0000 


3.85 + 1.30 


6+2 


3.0+0.5 


000 


680 


1.08 + 0.05 


22.4 


0.10 




0010 


2.63 + 0.06 


9+5 


1.4+0.4 


000 


684 


1.07 + 0.05 


22.3 


0.13 


0.9 + 0.1 


0000 


2.98 + 0.16 


17+3 


2.0+0.5 


000 


758 


1.20 + 0.05 


22.2 


0.26 


1.0 + 0.1 


0000 


3.06 + 0.11 


10+3 


2.2+0.7 


000 


764 


0.74 + 0.64 




0.13 




0000 


2.94 + 0.20 


9+3 


1.7+0.5 


000 


769 


1.03 + 0.02 


22.2 


0.17 




0000 


3.05 + 0.21 


11+4 


2.1+0.6 


000 


783 


1.18 + 0.08 


22.3 


0.17 




0010 


2.80 + 0.23 


15+3 


2.8+0.5 


000 


797 


1.07 + 0.05 


22.5 


0.17 


1.3 + 0.2 


0010 


2.89 + 0.28 


11+3 


1.6+0.5 


000 


798 


1.13 + 0.07 


22.2 


0.30 


0.8 + 0.1 


0000 


2.87 + 0.11 


10+3 


2.0+0.8 


000 


806 


1.05 + 0.05 


22.2 


0.21 


1.1+0.1 


0000 


2.99 + 0.19 


8+2 


1.8+0.5 


000 


820 


1.13 + 0.09 


22.3 


0.15 


0.7 + 0.1 


0011 


2.99 + 0.23 


6+1 


3.6+0.7 


000 


822 


1.10 + 0.06 


22.1 


0.16 


0.8 + 0.1 


0000 


3.02 + 0.08 


13+4 


1.6+0.5 


000 


823 


1.18 + 0.04 


22.0 


0.26 


0.9 + 0.1 


0000 


2.84 + 0.12 


12+2 


2.1+0.7 


000 


824 


1.19 + 0.05 


21.6 


0.24 


1.1+0.2 


0000 


2.77 + 0.08 


11+2 


4.0+1.1 


000 


826 


1.17 + 0.05 


22.0 


0.17 




0000 


2.58 + 0.32 


14±4 


3.0+0.9 


000 


828 


1.18 + 0.05 


22.6 


0.16 


1.0 + 0.1 


0010 


2.78 + 0.07 


15+4 


1.3+0.4 


000 


829 


1.25 + 0.08 


22.0 


0.17 


1.4 + 0.2 


0010 


4.00 + 0.69 


13+3 


2.2+0.3 


000 



Article number, page 22 of 27 



Husemann et al. : The CALIFA survey IL First public data release 



Table 6. continued. 



ID 


airmass 




Av 


seeing 


flags(R) 


SA 


S/N(/?5o) 


her 


flags(QC) 


832 


1.03 + 0.02 


22.1 


0.17 




0000 


2.92 + 0.27 


18+5 


2.1+0.6 


000 


833 


1.07 + 0.03 


22.2 


0.10 




0000 


2.99 + 0.69 


5+2 


3.4+1.2 


000 


835 


1.04 + 0.04 


22.3 


0.16 


1.0 + 0.1 


0000 


2.66 + 0.11 


20+5 


1.4+0.5 


000 


837 


1.16 + 0.05 


22.2 


0.14 


0.8 + 0.1 


0010 


2.91 + 0.13 


25+6 


1.6+0.5 


000 


840 


1.26 + 0.11 


22.1 


0.10 




0000 


3.71 + 1.47 


12+3 


2.1+0.9 


000 


843 


1.12 + 0.06 


22.4 


0.26 




0000 


2.74 + 0.15 


10+4 


1.9+0.4 


000 


845 


1.19 + 0.17 


21.7 


0.28 


0.8 + 0.1 


0000 


2.92 + 0.10 


11+4 


2.3+1.0 


000 


846 


1.18 + 0.08 


22.0 


0.25 




0000 


2.92 + 0.14 


9+2 


2.0+0.9 


000 


847 


1.03 + 0.03 


21.9 


0.32 


0.8 + 0.1 


1000 


2.95 + 0.27 


11+4 


2.3+0.8 


000 


848 


0.83 + 0.73 


14.7 


0.23 


1.0 + 0.1 


1000 


2.66 + 0.11 


16+3 


2.3+0.6 


000 


850 


1.05 + 0.02 


22.0 


0.19 


0.8 + 0.1 


0000 


3.00 + 0.10 


13+3 


2.2+0.5 


000 


851 


1.13 + 0.04 


22.1 


0.28 


0.9 + 0.1 


0000 


3.00 + 0.74 


8+2 


2.0+0.5 


000 


852 


1.40 + 0.10 


21.6 


0.19 


1.4 + 0.2 


0010 


3.34 + 0.39 


6+2 


2.1+0.6 


000 


854 


1.22 + 0.06 


22.0 


0.37 


0.9 + 0.1 


0001 


2.69 + 0.16 


10+4 


2.2+0.8 


000 


856 


1.06 + 0.04 


22.7 


0.24 


0.9 + 0.1 


1100 


2.69 + 0.22 


7+2 


3.2+0.8 


000 


857 


1.10 + 0.02 


22.4 


0.22 


0.9 + 0.1 


0000 


2.83 + 0.13 


16+4 


1.9+1.1 


000 


858 


1.28 + 0.17 


22.1 


0.23 


0.9 + 0.1 


0000 


2.97 + 0.27 


13+5 


2.4+0.8 


000 


859 


1.18 + 0.05 


22.2 


0.32 


0.7 + 0.1 


0000 


2.65 + 0.10 


14+3 


2.6+0.6 


000 


860 


1.05 + 0.03 


22.3 


0.27 


0.9 + 0.3 


0000 


2.78 + 0.31 


23+3 


2.4+0.6 


000 


863 


1.20 + 0.06 


22.3 


0.19 


0.9 + 0.1 


0000 


2.68 + 0.07 


13+4 


2.2+0.8 


000 


864 


1.05 + 0.02 


22.4 


0.20 




0000 


3.56 + 0.32 


14+3 


1.7+0.5 


000 


865 


1.02 + 0.02 


22.3 


0.26 


0.9 + 0.1 


0000 


2.94 + 0.14 


17+4 


1.7+0.6 


000 


866 


1.08 + 0.05 


22.1 


0.34 




0000 


2.95 + 0.06 


6+2 


1.7+0.4 


000 


867 


1.13 + 0.01 


22.1 


0.38 




0000 


3.61 + 2.32 


17+2 


2.4+0.6 


000 


872 


1.37 + 0.06 


21.7 


0.22 


1.0 + 0.1 


0000 


2.66 + 0.14 


7+2 


2.8+0.8 


000 


874 


1.13 + 0.05 


21.8 


0.30 


0.8 + 0.1 


0000 


3.29 + 0.47 


11+3 


2.7+0.4 


000 


877 


1.10 + 0.06 


21.7 


0.22 


1.0 + 0.1 


0000 


3.01 + 0.22 


8+3 


3.5+1.0 


000 


878 


1.36 + 0.06 


22.1 


0.18 


1.1 +0.1 


0010 


2.86 + 0.11 


5+7 


2.1+0.7 


000 


881 


1.17 + 0.12 


22.0 


0.30 


0.9 + 0.2 


0000 


3.04 + 0.19 


16+4 


2.0+0.7 


000 


883 


1.12 + 0.04 


22.1 


0.31 


0.9 + 0.1 


0000 


2.71 + 0.19 


3+1 


4.6+1.7 


000 


887 


1.05 + 0.01 


22.2 


0.29 


0.8 + 0.1 


0100 


3.03 + 0.42 


13+4 


2.8+1.0 


000 


888 


1.08 + 0.06 


21.9 


0.41 




0000 


2.87 + 0.14 


8+3 


2.5+0.8 


000 


890 


1.02 + 0.01 


21.8 


0.23 


0.9 + 0.1 


0000 


3.03 + 0.22 


13+3 


2.5+0.6 


000 


893 


1 on 1 r\ o 1 

1.37 + 0.21 


21.8 


0.15 


1.0 + 0.1 


0010 


2.8 / + U.14 


8+3 


2.4+0.4 


000 


896 


1.04 + 0.03 


22.3 


0.14 


1.4 + 0.3 


0110 


2.80 + 0.30 


11+3 


2.5+0.6 


000 


901 


1.20 + 0.11 


22.1 


0.25 


1.0 + 0.2 


0100 


3.09 + 0.15 


11+6 


3.5+1.7 


000 


902 


1.34 + 0.14 


21.7 


0.21 


1.1+0.1 


0000 


3.04 + 0.15 


18+3 


1.9+0.5 


000 


904 


1.29 + 0.09 


21.7 


0.22 


1.3 + 0.3 


0001 


3.03 + 0.25 


14+4 


2.5+0.8 


000 


935 


1.10 + 0.07 


22.3 


0.18 




0110 


2.43 + 0.04 


6+2 


2.5+0.7 


000 


939 


1.35 + 0.08 


21.9 


0.28 


0.9 + 0.1 


0100 


3.23 + 0.28 


9+3 


1.9+0.7 


000 



Notes. The meaning of the columns is the same as for Table [5] except that the sky surface brightness is measured in the B band and the QC flag 
(FLAG_FLUXCAL) is not defined for the V1200 data. 
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Appendix A: CALIFA data reduction pipeline 
improvements 

Appendix A.1: Detection of cosmics in single exposures 

Although the individual exposure times for the CALIFA data are 
relatively short, numerous pixels on the CCD may still be hit by 
cosmic rays. These pixels must be identified and properly han- 
dled to avoid artifacts in the final data products. In pipeline VI .2, 
cosmic rays in single exposures were detected using a simplified 
Laplacian edge detection scheme based on van Dokkum (2001 ). 
However, the performance was limited, with many cosmics re- 
maining undetected. For this reason we developed a new detec- 
tion algorithm (PyCosmic, Husema nn et al.|20T2l ) that achieves 
a much higher detection efficiency by combining a Laplacian 
edge detection scheme with a point spread function convolution 
approach. PyCosmic is applied to each raw frame just after it 
has been bias-subtracted and converted to electron counts. All 
detected cosmic rays are replaced by the median value of their 
surrounding pixels in PyCosmic. Nevertheless, we use the asso- 
ciated bad pixel mask to properly mask those pixels throughout 
the VI. 3 c pipeline as suggested in [Husemann et al.|p012[ ). 



NGC6478_p3 flexure shifts (x200) 



Appendix A. 2: Estimation and correction of instrument 
flexures 



The PMAS instrument suff'ers from significant flexures (" Rothl 
|et al.||2003] ), which introduces small shifts of the signal on the 
CCD in the dispersion and cross-dispersion direction depend- 
ing on the hour angle of the telescope. Pipeline VI. 2 invoked 
a pre-defined position list of ~30 prominent night-sky emission 
line across the CCD to estimate the absolute shifts in both di- 
rections. Subsequently, each raw image was resampled to a 
common reference frame if the off'set exceeded 0.1 pixel. This 
method was unstable because the absolute flexure ofl'sets can ex- 
ceed several pixels in which case the sky line positions might be 
miss-identified with lines in the neighboring fibers. 

In pipeline VI. 3c, we only estimate relative flexure ofl^sets 
with respect to the continuum and arc-lamp calibration frame 
of each object in a two- stage process. Relative off'sets in the 
cross-dispersion direction are measured by comparing the peak 
positions of fibers in each science frame close to the brightest 
night-sky emission lines with respect to their position in contin- 
uum lamp frame. After extraction of spectra, the wavelength of 
prominent sky lines is measured assuming the wavelength so- 
lution imposed by the arc-lamp calibration, which is compared 
with the expected wavelength of those lines. Any systematic 
deviation is attributed to flexure efl'ects and the wavelength solu- 
tion is corrected specifically for each individual science frame. 
These relative flexure ofl'sets are often small (< 0.5 pixel) be- 
cause of the relatively small changes in telescope position be- 
tween the science and calibration exposures for CALIFA. An 
example of the flexure pattern is shown in Fig. |A.l[ which indi- 
cates that the flexure shifts are nearly constant across the field. 
A global shift in dispersion and cross-dispersion is adopted for 
each frame. No resampling of the data is applied to correct for 
flexures in pipeline VI. 3c, because only the calibration informa- 
tion is adjusted. 

Appendix A.3: Subtraction of straylight 

An undesired feature of PMAS is an almost homogeneous back- 
ground signal across the CCD image. Scattered light within the 
spectrograph, highly variable dark currents, or issues with the 
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Fig. A.l. Example of the measured flexure shifts for the third V500 
pointing of NGC 6478. The black arrows indicate the shift in x and 
y direction in CCD detector coordinates for positions across the CCD 
in the science exposures with respect to the continuum and arc lamp 
calibration frames. The length of the arrows represent the shift in pixels 
but multiplied by a factor of 200 to be visible. The median shift of all 
measurements is shown by the red arrow. 



CCD controller could be potential origins and we simply refer 
to this feature as "straylight". We find that the straylight con- 
tribution to some fibers can exceed 15% in counts depending 
on their throughput and input signal. The strength of the stray- 
light aff'ects the S/N of an object because it increases the Poisson 
noise. Straylight was implicitly corrected in pipeline VI. 2 by 
the sky subtraction, because the sky spectra inherit almost the 
same straylight signal as all other spectra. This approach, how- 
ever, takes neither the subtle inhomogeneity of the straylight into 
account, nor assigns the correct error budget to the individual 
constituents, i.e. straylight, sky and object signal. 

In pipeline VI. 3c, we reconstruct the straylight signal using 
the few gaps between the fiber traces. First, the data is smoothed 
in the dispersion direction using a 5 pixel wide median filter to 
suppress the noise of individual CCD pixels. In a second step, we 
fit a sixth order polynomial to each cross-dispersion slice of the 
CCD, considering only pixels more than 6 pixels away from any 
fiber peak. An example is shown in Fig. | A. 2 [ (upper panel). Fi- 
nally, the resulting straylight map is smoothed by a 2D Gaussian 
kernel with a dispersion of 20 pixels to suppress any artificial 
high-frequency structure. Ex ampl es of the smoothed straylight 
distribution are shown in Fig. | A. 2 [ (lower panels) for the VI 200 
and V500 setup. This distribution is subtracted from the calibra- 
tion and science exposures before the fiber spectra are extracted. 
The straylight signal is of the order of 5-30 counts for the major- 
ity of the data, but in very few cases it reached several hundred 
counts or showed a complex pattern potentially due to problems 
with the cooling of the CCD. The computed Poisson noise of the 
images is unchanged so that the eff'ect of straylight is properly 
included in the error budget of the data. 
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Fig. A.2. Example of the straylight subtraction process. Top panel: 
Pixel counts along a CCD cross-dispersion cut of the science observa- 
tion. A sixth order polynomial fit only through the gaps between the 
fiber bundles is shown as the red line. Bottom panel: Straylight maps 
for the V500 and VI 200 data, respectively, reconstructed from the poly- 
nomial fit for each cross-dispersion slice of the CCD and smoothed by 
a 2D Gaussian filter with a 20 pixel width. 



Appendix A.4: Extraction of the fiber spectra 

The amount of cross-talk between the PPak fibers was esti- 
mated in pipeline VI. 2 during spectral extraction by means of a 
modified Gaussian suppression technique ( |Sanchez|2QQ6| ). This 
scheme achieved an accuracy of less than 1% for the cross-talk if 
the width of each fiber profile is a-priori known along the disper- 
sion axis. A constant FWHIM of 2.5 pixels in cross-dispersion 
direction was assumed in pipeline VI. 2, but it changes signifi- 
cantly across the CCD. 

In pipeline VI. 3c, we now measure the widths of the fiber 
profiles averaged over blocks of 20 fibers at each 50th spectral 
element along the dispersion axis using a minimization with 
the integrated fiber fluxes and a common FWHIM per block as 
free parameters. The position of fibers is fixed based on the in- 
formation from the continuum lamp (including any flexure oflT- 
set). Subsequently, the measured fiber widths are interpolated 
by a 5th-order polynomial along the entire dispersion axis. An 
example of the estimated FWHIM along the dispersion axis is 
shown in Fig. |A.3| 

For the spectrum extraction in pipeline VI. 3c, we use the 
optimal extraction algorithm (Home 198 6 ) using a parallelized 
code that takes advantage of multi-core processor systems. The 
positions and widths of each Gaussian fiber profile are fixed to 
the values of the previous measurements. The advantage of opti- 
mal extraction is that errors can be analytically propagated based 
on the individual pixel errors. Bad pixels are also withdrawn dur- 
ing the extraction, despite a higher associated error. However, 
the extracted flux would be unreliable if all three of the most 
central pixels of a fiber are simultaneously flagged as bad. Thus, 
the corresponding spectral resolution element is flagged as bad 
in that case to avoid potential artifacts in the data. 




dispersion direction [pixel] 

Fig. A.3. Example of the measured fiber profile width in cross- 
dispersion along the dispersion direction for each fiber. The measure- 
ments are performed in blocks of 20 fibers and interpolated along the 
dispersion direction with a fifth order polynomial function. 

Appendix A.5: Wavelength calibration and honnogenization of 
the spectral resolution 

The wavelength solution, as well as changes of the spectral res- 
olution with wavelength are measured for each fiber based on 
the HeHgCd calibration lamp exposure taken for each CALIFA 
dataset. This information is used to resample the spectra to a lin- 
ear grid in wavelength and to homogenize the spectral resolution 
to a common value along the dispersion axis using an adaptive 
Gaussian convolution. This scheme is unchanged in the current 
version of the pipeline, but now the flexure oflTsets in the disper- 
sion direction are explicitly included in the wavelength solution 
to avoid a further resampling step. 

The main goal of the VI. 3c pipeline is to properly process 
the errors and bad pixel masks during the reduction step. The 
errors are analytically propagated during the Gaussian convolu- 
tion and a IVLonte Carlo approach is used to estimate the noise 
vector after the spline resampling of the spectra. Both process- 
ing steps introduce some correlation between neighboring pixels 
in the data and the noise. This is not a serious problem unless 
the data are binned in the dispersion direction. The assumption 
of uncorrelated errors would lead to a higher S/N compared to 
what will be measured from the real spectra. Bad pixels are 
completely masked out during both processes. Given the loss 
of information during the spline resampling at these pixels, we 
expand the bad pixel mask by 2 neighboring pixels on both sides 
along the dispersion axis. Finally, the error value of bad pixels 
is set to a high value (~ 10^^ counts) and the actual data are re- 
placed with a linear interpolation along the dispersion axis of the 
nearest unmasked pixels. 

Appendix A.6: Fiber transmission correction 

Applying a correction for the underlying straylight signal 
(Sect. |A.3| ) also improves the quality of the fiber flat-fielding. 
This is particularly important for the vignetted areas at the cor- 
ners of the CCD (see S 12 for the details of the origin and eflfect 
of the vignetting), where the contribution of the straylight starts 
to dominate over the actual fiber signal. Consequently, we were 
able to reduce the threshold in the relative fiber transmission at 
which the data suflTer from very low S/N from 75% to 30% in 
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Fig. A.4. Relative sensitivity curves normalized at 6000 A from sim- 
ulated (noise free) standard star observations with PPak at a seeing of 
1.2'\ Single pointing observations were simulated for a grid of point- 
ing positions in a and 6 with respect to the standard star. Simulations 
were performed for airmasses of 1.5 (gray), 1.3 (green), 1.1 (blue), and 
1.0 (red) to take into account the change in the differential atmospheric 
refraction. The blue-to-red photometric calibration accuracy achievable 
is not sufficient even at low airmasses simply because of wavelength 
dependent aperture losses. 
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Fig. A.5. The relative blue-to-red ratio between SDSS spectra and ob- 
served CALIFA V500 galaxy spectra (matched to the SDSS fiber aper- 
ture) along wavelength. All CALIFA spectra shown were obtained dur- 
ing the photometric night (2012 January 20th) and flux calibrated using 
the single pointing standard star observation of that night. A fourth or- 
der polynomial functions is sufficient to describe the trend presumably 
originating from the aperture losses of the standard star observation. 



pipeline VI. 3c. Pixels below this threshold are flagged as bad 
pixels and set to zero value. This decreases the masked area 
caused by vignetting at the comers of the CCD. It also increases 
the reconstructed image quality because three low transmission 
fibers (55-70% of transmission) within the PPak FoV now re- 
main above this new transmission threshold and do not need to 
be masked out anymore. 

Appendix A. 7: Flux calibration 

The large fiber size and low filling factor of PPak can complicate 
obtaining an accurate sensitivity function for the whole wave- 
length range. The relative blue-to-red sensitivity function de- 
pends strongly on the actual pointing of the PPak fiber bundle 
with respect to the position of the standard star. The strength of 
this eflTect can be severe as inferred from a large suite of simu- 
lated standard star observations (see Fig. |A.4| ). The reason for 
this is that the light of a standard star covered by a PPak fiber 
inevitably changes with wavelength because of atmospheric dis- 
persion, thor oughly dis cussed in|Garcia-Benito| ( |20Q9| ), [Garcia- 1 
Benito et alT] ( i201Q) and |Rosales-Qrtega|pQ09| ). These classical 
aperture losses occur even if the standard star has been centered 
on a single fiber as accurately as possible for panchromatic ob- 
servation as performed with the V500 setup. 

We therefore empirically calibrated these aperture losses by 
comparing the CALIFA spectra with available SDSS spectra. 
We selected two photometric nights (shortly before and after 
mirror re-coating in September 201 1) which had the best condi- 
tions: (i) the standard star was almost centered on a single fiber, 
(ii) the seeing conditions were good (<1.5''), and (iii) the airmass 
was sufl&ciently low (<1.3). After extracting the calibration star 
spectra we corrected them for the atmospheric extinction along 
wavelength (San chez et al.||2QQ7a| ) considering the airmass and 
the monitored V band extinction at the time of observation. The 
corresponding sensitivity function was then used for flux cali- 
bration of the science data. In contrast to pipeline VI. 2, we do 



not smooth the sensitivity function with a Gaussian kernel, but 
model it with a high-order polynomial to better preserve its shape 
while clipping outliers and residuals caused by stellar absorption 
line mismatches. 

SDSS DR7 spectra are available for -40% of the CALIFA 
galaxies observed during those two nights. We extracted the 
spectrum within a radius of 1.5'' around the galaxy center tak- 
ing fractional pixel coverage into account to properly compare 
the spectrophotometry of the SDSS spectra and CALIFA data. 
The ratio of the 3'' -d iame ter aperture SDSS and CALIFA spec- 
tra are shown in Fig. |A.5| for one of the nights. The significant 
red-to-blue deviation along the wavelength is attributed purely to 
aperture losses of standard star observations as discussed before. 
We corrected the initial sensitivity functions of the two nights 
correspondingly and used the mean of both as a master sensitiv- 
ity function for the calibration of all the CALIFA data distributed 
in DRl. The quality of this approach is evaluated in Sect. [5.4.2 



Appendix A.8: Sky subtraction 

A high S/N sky spectrum is obtained by combining the 36 sky 
fibers of PPak. In pipeline VI. 2, we rejected outlier sky spec- 
tra using a 2cr clipping. This approach turned out to be unstable 
when a bright field star or neighboring galaxy filled an entire sky 
fiber bundle consisting of 6 fibers. Thus, we decided to simply 
take the mean of the 30 faintest sky fibers for pipeline VI. 3c, 
which eflTectively avoids this problem except if more than one 
sky fiber bundle would be aflfected. In such a rare case it would 
be difliicult anyway to obtain a reliable sky spectrum from the 
few remaining sky fibers. Another advantage of this procedure 
is that the error vector for the mean sky spectrum is well-defined 
and easily propagated through the sky subtraction step. One 
caveat is that uncertainties in the wavelength solution and the 
spectral resolution matching are not included in the uncertainty 
of the sky subtraction. 
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Appendix A.9: Image reconstruction and DAR correction 

Since pipeline VI. 2 a flux-conserving inverse-distance weight- 
ing scheme (see S 12 for details) is used to reconstruct a spatial 
image from the three dithered PPak pointings with a sampling of 
V\ The efl'ect of diff'erential atmospheric refraction (DAR) on 
the data can be empirically measured from the spatially resam- 
pled datacube by tracing the centroid of the galaxy on the image 
plane along the wavelength axis. Common practice is often to 
resample each monochromatic image slice of the datacube again 
to align the galaxy centroid to a common reference position. 

It is crucial in pipeline VI. 3 to avoid this additional resam- 
pling step for a clean error propagation. Therefore, we adopt 
a two-stage iteration, (i) reconstruct the datacube and estimate 
the DAR eff'ect, (ii) reconstruct the datacube again where the 
position of the fiber against the regular grid is changed accord- 
ing to the DAR off'set empirically measured in the previous step. 
Although the error vector is well-defined for each spectrum in 
the reconstructed cubes, the inverse-distance weighting scheme 
introduces a strong correlation of the signal in the spatial dimen- 
sion which is discussed in Sect. l3.3l in more detail. 

Bad pixels are completely masked out during the weighting 
scheme so that in principle a full image can be reconstructed. 
However, a pixel value cannot be appropriately reconstructed if 
the contrast between the expected flux and the measured flux in 
the neighboring fibers is too high. We therefore flag each pixel 
in a reconstructed image as bad if the relative distance-based 
weight exceeds 5% for a fiber masked as bad at a given spec- 
tral resolution element. This value was chosen by experience to 
properly flag apparent artifacts along the wavelength axis caused 
by this efl'ect. 
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